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FOREWORD 


At  the  present  time,  thousands  and  thousands  of  people  around  the  world  deal 
with  ice,  snow  and  permafrost.  They  are  scientists,  educators,  engineers,  navigators, 
meteorologists  and  others.  While  a  small  fraction  of  these  people  contribute  to  the 
knowledge  base  in  ice  physics,  all  of  them  use  knowledge  from  it  frequently. 
Moreover,  successful  applied  research  is  based  upon  fundamental  science — one 
more  reason  for  ice  specialists  to  have  a  textbook  on  ice  physics  on  their  desks. 

The  first  modem  ice  physics  text  was  Fletcher's  book  on  The  Chemical  Physics  of  Ice 
(1970).  Fletcher's  book  is  in  typical  textbook  format:  it  is  reasonably  brief  and  easy 
to  understand.  He  touched  on  a  few  of  the  most  important  topics,  but  not  all  of  them. 

The  most  comprehensive  book  on  ice  physics  to  date  was  published  by  Hobbs  in 
1974.  Hobbs  considered  almost  all  of  the  basic  aspects  of  ice  as  understood  at  that 
time.  Moreover,  he  described  and  compared  several  (sometimes  opposing)  view¬ 
points.  This  fundamental  and  rather  large  (837  pages)  book  is  commonly  known  as 
the  "Ice  Bible"  by  specialists  in  the  field.  In  1974  and  1975,  two  CRREL  Monographs 
on  ice  were  produced  by  John  Glen.  These  were  briefly  and  clearly  written  and 
reviewed  almost  all  ice-related  subjects.  This  work  was  (and  in  some  respects  still  is) 
a  magnificent  introduction  to  ice. 

Finally,  in  1981  Maeno  wrote  a  simple  popular  book  for  the  express  purpose  of 
attracting  people's  attention  to  the  subject. 

During  the  past  20  years,  a  significant  amount  of  new  experimental  and  theoreti¬ 
cal  work  has  appeared,  dramatically  changing  our  views  on  ice  physics.  As  a  result, 
we  are  now  able  to  formulate  physical  laws  using  more  simple  and  direct  methods. 
We  have  found  some  of  the  physical  models  used  in  the  past  to  be  completely  wrong. 
The  physics  of  ice  is  a  much  better  developed  subject  than  it  was  20  years  ago. 

For  the  above  reasons,  we  feel  the  time  is  ripe  for  a  contemporary  book  on  ice 
physics,  incorporating  the  known  and  proven  with  almost  20  years  worth  of 
material  not  covered  by  previous  works. 

We  have  tried  to  prepare  a  "readable"  book,  and  not  one  that  requires  the  reader 
to  be  a  uniquely  educated  person.  It  is  our  intent  to  present  the  material  in  such  a  way 
that  any  reader  attracted  by  the  title  Ice  Physics  will  be  able  to  comprehend  it.  This 
is  quite  difficult  for  a  book  dedicated,  not  to  a  particular  field  of  knowledge,  but  to 
a  specific  material.  Indeed,  for  ice  it  means  we  have  to  consider  a  wide  variety  of 
subjects,  including  quantum  chemistry,  solid  state  physics,  the  theory  of  elasticity, 
ionic  conductivity,  synchrotron  x-ray  topography,  crystal  growth,  the  physics  of 
surfaces  and  more. 

The  primary  goal  is  to  produce  as  simple  a  book  as  possible  without  sacrificing 
scientific  accuracy.  Experimental  facts,  physical  ideas  and  theories  will  be  strongly 
organized  and  bound  together  cohesively.  The  reader  will  be  introduced  to  a  wide 
variety  of  material  on  a  step  by  step  basis.  Then  the  picture  will  be  whole. 

To  accelerate  materials  publication,  this  book  will  appear  first  in  the  form  of  a 
series  of  joint  CRREL-Dartmouth  reports,  later  to  be  published  in  CRREL's  Mono¬ 
graph  series,  on: 

1.  The  structure  of  ordinary  ice 

Part  I:  "Ideal"  structure  of  ice.  Ice  crystal  lattice 
Part  II:  Defects  in  Ice 

Volume  1:  Point  defects 

Volume  2:  Dislocations  and  planar  defects 

2.  Electrical  properties  of  ice 

Part  I:  Conductivity  and  dielectric  permittivity  of  ice 
Part  II:  Advanced  topics  and  new  physical  phenomena 

3.  Optical  properties 

4.  Electro-optical  effects  in  ice 


iii 


5.  Theimal  properties 

6.  Mechanical  properties  of  ice.  Elasticity  and  anelastic  relaxation.  Plastic 
properties.  Fracture  of  ice 

7.  Electromechanical  effects  in  ice 

8.  Surface  of  ice 

9.  Other  forms  of  ice  and  their  properties 

10.  Ice  in  space 

11.  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  convenient  to  the  authors.  The  present 
report  is  the  third  in  the  series. 
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Structure  of  Ordinary  Ice  Ih 
Part  II:  Defects  in  Ice 
Volume  1:  Point  Defects 


VICTOR  F.  PETRENKO  AND  ROBERT  W.  WHITWORTH 


INTRODUCTION 

Unlike  most  previous  authors  of  books  on  ice 
physics  (Hobbs  1974,  Fletcher  1970,  Glen  1974),  we 
decided  to  isolate  all  kinds  of  defects  in  ice  into 
separate  reports.  Such  a  division  reflects,  first,  the 
important  role  that  defects  play  in  materials  phys¬ 
ics.  While  most  of  the  physical  properties  of  solids 
are  structurally  sensitive  ones,  in  ice  some  of  them 
are  entirely  determined  by  specific  defects.  For 
example,  ion  and  Bjerrum  defect  fluxes  are  re¬ 
sponsible  for  electrical  properties  and  inelastic 
relaxation  in  ice  (see,  e.g.,  Petrenko  and  Ryzhkin 
1984c).  Second,  it  is  in  the  study  of  defects  (dislo¬ 
cations,  point  defects,  grain  boundaries)  that  sig¬ 
nificant  progress  has  been  achieved  since  the  last 
monograph  on  ice  physics  was  published  (Hobbs 
1974).  Because  of  this,  special  attention  will  be 
given  to  the  results  of  studies  from  the  last  two 
decades. 

Following  a  tradition  in  materials  science,  we 
divide  consideration  of  defects  into  three  parts: 
point  defects,  dislocations  and  planar  defects.  In 
this  report  we  discuss  the  structure  and  physical 
properties  of  zero-dimensional  defects,  known  as 
point  defects.  In  the  second  report,  entitled  Struc¬ 
ture  of  Ordinary  Ice  Ih;  Part  II:  Defects  in  Ice;  Volume 
2:  Linear  and  Planar  Defects,  we  consider  disloca¬ 
tions  that  are  one-dimensional  defects,  and,  fi¬ 
nally,  grain  boundaries  and  stacking  faults  (two- 
dimensional  defects)  will  be  described.  The  struc¬ 
ture  and  properties  of  another  two-dimensional 
"defect"— the  surface — are  treated  in  a  separate 
report. 

GENERAL  CONSIDERATION 
OF  POINT  DEFECTS 

We  begin  our  review  of  defects  in  ice  with  the 
smallest  and  simplest  ones.  These  are  the  point 
defects  that  are  usually  vacancies,  interstitials, 
impurity  atoms,  bivacancies,  Schottky  and  Frenkel 
pairs,  etc.,  but  in  ice  also  include  special  kinds  of 


defects  associated  with  the  disorder  of  the  hydro¬ 
gen  atoms.  The  term  "point"  does  not,  of  course, 
mean  that  the  defects  do  not  have  a  size,  but  that 
they  are  located  at  a  specific  point  within  the 
crystal.  They  are  usually  located  within  one  or 
two  unit  cells,  although  elastic  deformation  pro¬ 
duced  by  the  defects  can  spread  to  much  larger 
distances. 

Point  defects  in  ice  (as  in  other  solids)  play  a 
fundamental  role  in  many  physical  processes  such 
as  diffusion,  anelastic  relaxation,  electrical  con¬ 
ductivity,  dielectric  permittivity  and  optical  ab¬ 
sorption.  They  differ  from  one  another  in  struc¬ 
ture,  concentration  and  mobility,  but  in  this  intro¬ 
duction  we  will  consider  some  general  results 
applicable  to  them  all.  These  are  the  relationships 
between  concentration,  formation  energy,  jump 
rates,  diffusion  coefficient  and  mobility. 

Since,  for  a  given  number  of  molecules  the 
arrangement  with  minimum  energy  is  perfect, 
with  no  defects,  some  additional  energy  Ef  is  re¬ 
quired  to  produce  a  defect.  For  the  majority  of 
point  defects  in  ice,  Ef  is  of  the  order  of  0.5  eV. 
These  values  are  much  greater  than  the  aver¬ 
age  energy  of  thermal  vibration,  fcgT  =  0.023  eV 
(-10°C),  so  that  at  first  sight  the  formation  of  these 
defects  is  unjustified.  However,  according  to  ther¬ 
modynamics  at  constant  temperature  and  con¬ 
stant  volume,  thermal  equilibrium  is  attained  at 
minimum  free  energy  F  rather  than  at  minimum 
internal  energy  U,  where 

F  =  U-TS.  (1) 

The  appearance  of  defects  increases  the  internal 
energy  U,  but,  because  of  the  large  number  of 
different  ways  in  which  the  defects  can  be  located 
in  the  crystal,  their  entropy  increases  more  rap¬ 
idly,  and  the  minimum  of  F  occurs  at  a  finite  defect 
concentration. 

We  can  illustrate  this  point  with  a  simple  ex¬ 
ample.  Consider  a  crystal  lattice  consisting  of  N 
identical  points  (atoms,  molecules  or  elementary 
cells).  A  point  defect  with  energy  of  formation  Ef 


could  appear  with  equal  probability  at  any  of  the 
points.  Then,  the  free  energy  attributable  to  the 
formation  of  n  defects  will  be 

Fd  =  Ef/J-TSc  (2) 

where  Sc  is  the  configurational  entropy  of  the  n 
defects  in  the  lattice 


Sc  =  EglnW 


(3) 


and  W is  the  number  of  arrangements  of  n  identical 
defects  over  these  N  identical  sites 


W  = - m - . 

n\(N  -  n) ! 
At  equilibrium 

^  =  0 


(4) 

(5) 


and  substituting  eq  2-4  into  eq  5  and  using  the 
Stirling  formula — ln(N!) = N  In  N-N  for  N  » 1 — we 
obtain 


/i 

N-n 


--  exp 


(6) 


activated  hops  over  these  barriers.  Let  us  consider 
the  motion  of  some  electrically  charged  defect 
with  charge  q,  which  we  will  refer  to,  for  brevity,  as 
an  "ion." 

The  simple  case  of  the  one-dimensional  motion 
of  such  an  ion  is  shown  in  Figure  1.  To  move  one 
interatomic  distance  in  the  direction  of  the  field  E, 
the  positive  ion  would  need  an  activation  energy 
Um-qEa/2t  and  to  move  against  the  field,  it  would 
need  Um  +  qEa/2.  The  frequency  of  ion  jumps  in 
the  direction  of  the  field  is  expressed  by 


f\ 


I  Um  -  qaE/2 1 
v  exp  [ - ) 

I  *b  T  I 


(9) 


and  against  the  field  by 


f2  =  v  exp 


I 


Um  -  qaE/2  | 

k^T  j 


(10) 


where  v  is  related  to  the  frequency  of  oscillations 
of  the  ion  in  its  potential  well.  Since  every  jump 
moves  the  ion  over  the  distance  a,  the  mean  drift 
velocity  of  the  ion  in  the  direction  of  the  field 
becomes 


Assuming  Ef »  k^T,  we  obtain 

H  =  JVexp|-pLj-  W 

To  be  more  precise,  the  appearance  of  a  point 
defect  not  only  changes  the  configurational  en¬ 
tropy,  it  also  affects  the  vibrational  entropy  be¬ 
cause  of  a  slight  change  in  the  elastic  lattice  vibra¬ 
tion  spectrum.  This  vibrational  entropy  change  Sf 
would  lead  to  an  additional  factor  of  exp(SfAB)  hi 
eq  7  so  that 


_JL). 

(8) 

UB/  l 

kf?) 

Calculations  of  S f for  defects  are  quite  difficult, 
and  are  far  beyond  this  approach.  We  indicate  here 
only  that,  according  to  experimental  results,  the 
ratio  SfAc  may  not  be  small  and  must  be  taken  into 
account. 

The  majority  of  point  defects  in  ice  have  a  high 
mobility  at  "not  too  low"  (T>  100  K)  temperatures, 
migrating  from  one  lattice  site  to  another  by  jump¬ 
ing.  The  equilibrium  positions  of  point  defects  in 
the  lattice  are  separated  by  potential  energy  barri¬ 
ers  of  height  Um,  and  motion  occurs  by  thermally 


<Ud>  =  a  {f\  -f2)  =  2v  a  exp 


Urn) 

AbT| 


sinh 


1  qEa 
,2  kBT 
(11) 


Since  usually 

q  E  a  «  2  k^T  (12) 

we  can  substitute  sinh  with  its  argument  in  eq  11 


<ud> 


f£) 


exp 


(13) 


Thus,  the  ion  mobility,  which  is  defined  as  the  ratio 
between  drift  velocity  and  electric  field  strength  E, 
will  be 


b  = 


(14) 


and  is  independent  of  E.  Um  is  called  the  activation 
energy  of  migration. 

The  quantity 


1 

vh  =  -  v  exp 


UJl 

As  TJ 


(15) 
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Figure  1.  Potential  energy  of  an  ion  in  a  crystal  lattice  W  as  a 
function  of  coordinate  x  with  (E  >  0)  and  without  (E  =  0)  an 
electric  field-,  a  is  an  interatomic  distance. 


represents  the  probability  of  jumping  or  the  mean 
number  of  jumps  in  a  given  direction  per  second, 
and  the  diffusion  coefficient  of  the  ions  is  D  =  Vh«2. 
For  chaotic  thermal  motion  of  this  kind,  the  mobil¬ 
ity  |i  and  the  diffusion  coefficient  D  are  always 
related  by  the  Einstein  relation 


so  that 

D  =  (va2)expj--^J.  (17) 

Equation  17,  which  does  not  contain  q,  is  of 
course  also  valid  for  neutral  defects  such  as  vacan¬ 
cies  and  interstitials  and  for  neutral  atoms  and 
molecules  dissolved  in  ice.  When  the  above  simple 
formulae  are  for  some  reason  not  able  to  describe 
the  motion  or  concentration  of  some  specific  de¬ 
fect  in  ice,  we  shall  perform  an  additional  analysis. 

MOLECULAR  DEFECTS 

The  principal  molecular  point  defects  in  ice  are 
vacancies  and  interstitials.  A  vacancy  is  a  site  from 
which  a  water  molecule  is  missing,  and  an  intersti¬ 
tial  is  an  extra  molecule  inserted  in  a  space  be¬ 


tween  the  normal  molecules  ofthestructure.  These 
defects  can  clearly  be  formed  in  pairs  by  the  dis¬ 
placement  of  a  molecule  off  its  normal  lattice  site 
to  an  interstitial  position,  but  they  can  also  be 
formed  independently.  For  example,  a  vacancy 
appears  when  we  remove  a  water  molecule  from 
its  regular  position  and  place  it  on  the  surface. 
Dislocations  and  grain  boundaries  can  also  emit  or 
absorb  vacancies  and  interstitials. 

In  thermal  equilibrium  at  a  given  temperature, 
the  defect  concentrations  are  determined  by  their 
formation  energy  Ef  (and  any  associated  entropy 
Sf).  To  a  first  approximation,  the  energy  required 
to  form  a  vacancy  is  equal  to  the  energy  required 
to  remove  a  molecule  from  the  bulk  of  the  crystal 
and  place  it  at  a  kink  on  a  surface  step.  In  so  doing, 
the  number  of  hydrogen  bonds  per  molecule  is 
reduced  from  four  to  two.*  Subsequently,  there 
will  be  an  elastic  relaxation  of  the  lattice  around 
the  vacant  site,  which  slightly  lowers  the  energy. 

The  formation  of  an  interstitial,  on  the  other 
hand,  involves  the  transfer  of  a  molecule  from  a 


’Devlin  (1992a)  recently  found  from  spectroscopic  data  that 
the  surface  molecules  of  crystalline  ice  are  three-coordinate 
bonded,  i.e.,  there  is  only  one,  not  two,  dangling  O-H  bond  per 
molecule.  Nevertheless,  if  a  crystal  is  built  up  by  adding 
molecules  to  the  surface,  the  number  of  bonds  is  two,  irrespec¬ 
tive  of  the  number  of  dangling  bonds  at  the  surface. 
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kink  site  on  the  surface  to  an  interstitial  position. 
This  again  breaks  t*ib  hydrogen  bonds  and  in¬ 
volves  a  short-range  fepulsive  interaction  among 
nearest  neighbors  ant  an  elastic  distortion  of  the 
surrounding  lattice.  These  considerations  suggest 
that  the  energy  of  formation  of  these  defects  will  be 
of  about  the  same  m%nitude  as  the  cohesive  en¬ 
ergy  per  molecule,  vVfcich  is  known  to  be  0.6  eV. 

More  accurate  calculations  using  computer 
simulations  were  mate  by  Cotteril  et  al.  (1973).  In 
their  work  they  used®  realistic  potential  of  inter¬ 
action  between  wafcar  molecules  introduced  by 
Ben-Naim  and  Stillinger  (1972).  Their  calculations 
were  only  restricted  fy  their  assumption  that  the 
molecules  were  absolutely  rigid  and  the  small  size 
of  the  molecular  clualer  used  in  their  model.  The 
energies  of  molecula^ioint  defects  (and  also  of  the 
Bjerrum  L-  and  D-difects  to  be  described  later) 
calculated  by  Cotteai  et  al.  are  summarized  in 
Table  1. 

Water  molecules  'Ian  diffuse  through  ice  by 
either  a  vacancy  or  afr  interstitial  process.  To  be 
specific,  let  us  consider  such  self-diffusion  as  oc¬ 
curring  by  the  motio#  of  vacancies.  The  diffusion 
flux  of  vacancies  yv  rfgiven  by  the  formula 

jv  =  ~DV  grad((#j  (18) 

-I 

where  Dv  is  the  diffdkion  coefficient  of  vacancies 
defined  in  eq  17  andfcv  is  the  vacancy  concentra¬ 
tion.  A  particular  winter  molecule  in  the  ice  can 
migrate  by  vacancy*ldiffusion  only  if  one  of  its 
nearest-neighbor  ladlce  sites  happens  to  be  va¬ 
cant.  Consequently, 4he  coefficient  of  self-diffu¬ 
sion  of  molecules  is  much  smaller  than  the 
coefficient  of  vacanaj)  diffusion 


Table  1.  Defect  tarnation  energy  (eV)  in  ice  (after 
Cotteril  etal.  1973)." 

£ 

Cluster  size  and  clastic  relaxation 


Defect  16,  unrchdkd  16,  relaxed  43 2,  unrelaxed 

- 4 - 

Cohesive  0.466  4  —  0.465 

energy  , 

Vacancy  0.466‘^  0.447  0.465 

L-defect  0.467  €  0.388  0.403 

it 

D-defect  0.491  J  0.407  0.428 

Interstitial  0.616*  0.528  0.695 


D$  =  Dv  (19) 

«h20 

In  practice,  the  self-diffusion  coefficient  of  a  mate¬ 
rial  can  be  determined  using  a  tracer  technique, 
provided  that  the  material  possesses  suitable,  usu¬ 
ally  radioactive,  isotopes.  For  ice,  three  isotopes 
have  been  used  in  such  diffusion  measurements: 
deuterium,  2H;  tritium,  3H;  and  lsO.  The  general 
procedure  is  to  place  a  known  concentration  of  the 
isotope  onto  one  end  of  a  block  of  ice,  which  is  then 
left  at  a  fixed  temperature  for  a  certain  period  of 
time.  The  distribution  of  the  isotope  along  a  par¬ 
ticular  axis  is  then  measured  experimentally.  Table 
2  and  Figure  2  summarize  the  results  from  such 
experiments. 

If  we  substitute  eq  8  and  17  into  eq  19,  we  obtain 


Ds  -  exp  My 


in  which 


Oso={va^)exp 


£s  =  £f+Lim.  (22) 

Several  trends  in  the  experimental  data  in  Table  2 
are  of  interest.  First,  the  diffusion  coefficients  of 
oxygen  and  hydrogen  isotopes,  as  well  as  their 
activation  energies,  are  approximately  equal.  This 
is  major  evidence  that  in  ice  oxygen  and  hydrogen 
diffuse  cooperatively  by  the  diffusion  of  whole 
water  molecules.  The  scatter  in  Ds  and  £s  (rela¬ 
tively  small)  could  be  caused  by  different  degrees 
of  imperfection  of  the  ice  used  .For  the  maj  ority  of 
solids,  it  has  been  established  that  the  diffusion 
rate  along  grain  boundaries  or  along  dislocation 
lines  considerably  exceeds  the  bulk  rate.  Hence, 
one  could  expect  larger  values  of  Ds  in  polycrystal¬ 
line  ice  and  in  ice  with  high  dislocation  density 
(which  was  unknown  in  the  papers  cited). 

The  key  question  in  the  study  of  self-diffusion 
in  ice  is  about  whether  the  process  occurs  by  the 
motion  of  vacancies  or  interstitials.  Since,  in  the 
majority  of  crystals  vacancies  dominate  interstitials, 
self-diffusion  in  ice  was  for  a  long  time  also  ex¬ 
plained  in  terms  of  vacancy  diffusion.  However, 
experiments  by  Higashi,  Hondoh  and  others  dur¬ 
ing  last  decade  have  clearly  proven  that  the  domi¬ 
nant  process  in  diffusion  in  ice  is  that  ascribable  to 
interstitials.  Their  experiments  involved  the  care¬ 
ful  study  of  the  growth  and  annealing  of  disloca¬ 
tion  loops  produced  by  rapid  cooling. 
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Table  2.  Self-diffusion  coefficient  Ds  and 

Self  diffusion 


s  ctnu  its  activation  energy  E^in  ice. 

* 

* 


Measurement 

coefficient  us 
tm?  s~])  at 
-10°C 

Activation 
energy  Es  (eV) 

¥ 

Reference 

Oxygen1!)?  in  polycrystalline  ice. 

1.0  X  16-14 

— 

—fr — - 

Kuhn  Thiirkauf  (1958) 

Deuterium  in  polycrystallirie  ice. 

(at-i:0oCj 

i:6  X  It)  « 

___ 

t- 

Kuhn  Thurkauf  (1958) 

Tritium  perpendicular  to  c-Sxis  in  artificial  single 

(at-1.8°Cj 

i:5i  X  10-15 

6:58  ±  0.05 

* 

Denge^tnd  Riehi  (1^63) 

Tritium  perpendicular  and  parallel  to  c-axis  in 
natural  single  crystals.  Diffusion  time  of 

2:45  x  10-is 

0.6710.09 

Itagak|(1964) 

*• 

7  days. 

Tritium  parallel  to  c-axis  in  natural  single 

2.12  x  10  is 

0.54 1  0.09 

itagak>(1966, 1967) 

Tritium  perpendicular  to  c-axis  in  natural  single 

2.42  x  10  15 

0.65 1 0.05 

ltagakl(1966, 1967) 

crystals.  Diffusion  time  greater  than  2  days. 
Tritium  parallel  to  c-axis  in  artificial  single  crystals. 

2.15  X 10-15 

0.63 1. 0.03 

Biicksltai.  (1966) 

Tritium  perpendicular  to  c-aXis  in  artificial  single 

2.40  x  10  15 

0:63  ±  0:03 

Blicks  at  ai.  (1966) 

crystals. 

Oxygen-18  in  artificial  single  crystals. 

3.1  x  10-15 

0.68.10.13 

Deiibakas  et  ai.  (1966) 

Den  gel  et  ai.  (1966) 

Tritium  parallel  tti  c-axis  in  artificial  single  crystals. 

2.5  x  10-15 

0.63 1  0.05 

Tritium  parallel  to  c-axis  in  artificial  single  crystals. 

1.44  x  1(3-15 

0.62  ±  0.03 

Ramsidar  (1967) 

Tritium  perpendicular  to  c-axis  in  artificial  Single 

1.68  x  10-15 

0.62 1 0:07 

RamSeisr  (1967) 

Tritium  parallel  to  c-axis  in  natural  single  crystals. 

1.54  x  lO-D 

O.ei  ±  0.02 

Ramsiier  (i967) 

Tritiuin  perpendicular  to  c;axis  in  natural 
single  crystals. 

1:65  x  10-15 

0.62  ±  0.04 

Ramstdar  (1967) 

iS 

Consider  what  happens  in  the  crystal  during 
rapid  cooling.  The  term  "rapid"  means  that  the 
Coolihg  takes  place  during  a  period  much  less  than 
that  required  for  the  new  equilibrium  defect  con¬ 
centrations  to  be  established.  When  the  tempera¬ 
ture  is  lowered,  equal  numbers  of  vacancies  and 
interstitials  disappear  rapidly  owing  to  their  mu¬ 
tual  annihilation.  Then  some  fraction  of  defects 
that  are  predominant  must  ''die"  in  a  slower  and 
more  difficult  way — for  instance,  by  diffusing  to¬ 
wards  the  surface  or  concentrating  into  disk  clus¬ 
ters,  as  shown  in  Figure  3.  The  edges  of  these  disks 
are  edge  dislocation  loops,  with  the  Burgers  vector 
normal  to  the  plane  of  the  disk.  If  loops  of  a-type 
are  predominant,  then  the  major  defects  are  vacan¬ 
cies;  if  loops  bf  b-type  are  predominant,  then  the 
major  defects  are  interstitials. 


Dislocation  loops  (Fig.  can  expand  or  cbm 
tract  because  of  diffusion  df-point  defects  toward 
(or  from)  the  loops,  caused  by  a  change  in  temper¬ 
ature  Or  applied  pressure.  Monitoring  the  growth 
or  contraction  of  dislocation  loops  by  means  of  X- 
ray  topography  in  ice  monpcrystals,  Higashi  and 
Hondoh  with  co-authors  determined  the  diffusion 
coefficient  bf  interstitials  i»  ice  Dj,  the  formation 
energy  of  interstitials  Eas,  the  activation  energy  of 
their  motion  Ujm,  the  activation  volumes  of  inter¬ 
stitials  Vii  arid  of  their  migration  Vim,  and  absolute 
values  bf  the  interstitial  concentration  hi,  and  the 
entropy  SvAb-  Table  3  contains  ail  of  these  data. 

According  to  these  resalts,  the  activation  en¬ 
ergy  bf  self-diffusion  in  ice  $%  =  Eas  +  Ejm  =  0.56  eV, 
and  the  absolute  value  of  the  self-diffusion  coeffi¬ 
cient  is  irt  good  agreement  with  results  obtained 
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TEMPERATURE  (°C) 


Figure  2.  Comparison  of  the  self-diffusion  co¬ 
efficients  D  in  ice  obtained  by  the  tracer  tech¬ 
nique  (3H  and  lsO)  and  calculated  from  the 
diffusion  coefficients  of  self-interstitials  ob¬ 
tained  by  the  X-ray  topography  (solid  line) 
(after  Ogtiro  et  al.  19S8). 


103/T  (K-1) 


(a) 


<-b 


Figure  3.  Two  types  of  the  prismatic  dislocation  loop — 
(a)  vacancy  type  and  (b)  interstitial  type.  Schematic 
atomic  arrangement  at  the  top  corresponds  to  edge 
dislocation  loops  of  different  signs  (after  Oguro  et  al. 
1988). 
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Table  3.  Self-diffusion  interstitial  parameters  in  ice  at  0°C 
(after  Hondoh  1992).  Vmoi  is  a  molecular  volume  in  ice 
(VW  =  3.3  x  10-23  cm3). 

n,/N  Ejjs  (eV)  SdAb  D i(wj2/s)  Uj'm  (cV)  Vif/Vmol  Vim  /v mol 
2.8x10-6  0.40  4.9  2x10-9  0.16  -0.26  0.40 


earlier  by  the  radioactive  tracer  technique  (Fig.  2). 
We  refer  readers  interested  in  the  details  of  an 
elegant  experimental  technique,  employed  by 
Higashi  and  Hondoh  et  al.,  to  their  review  papers 
(Hondoh  et  al.  1987,  Oguro  et  al.  1988,  Hondoh 
1992),  where  one  can  find  references  to  their  origi¬ 
nal  papers. 

Thus,  from  measurements  of  self-diffusion  in 
ice  it  is  not  possible  to  determine  any  of  the  char¬ 
acteristics  of  vacancies,  since  the  predominant 
defects  are  interstitials  (at  least  at  T  >  -50°C).  The 
only  studies  in  which  activation  parameters  of 
vacancies  in  ice  were  possibly  determined  were 
those  in  which  a  positron  annihilation  technique 
was  used  (Mogensen  and  Eldrup  1978,  Eldrup  and 
Morgensen  1978).  When  ice  is  exposed  to  positrons, 
thermalized  positrons  and  positronium  atoms 
(electron  +  positron)  are  captured  mainly  by  va¬ 
cancies.  Therefore,  their  lifetimes  depend  on  the 
concentration  and  mobility  of  vacancies.  It  is  ex¬ 
tremely  difficult  to  estimate  the  absolute  accuracy 
of  this  method  since  it  is  based  on  a  considerable 
number  of  assumptions  and  approximations.  The 
vacancy  migration  energy  found  in  these  papers 
was  0.34  ±  0.07  eV  and  its  formation  energy  was 
found  to  be  between  0.2  and  0.35  eV. 


PROTONIC  DEFECTS 

As  explained  in  the  earlier  reports  on  the  Struc¬ 
ture  of  Ordinary  Ice  Ifo  Part  I:  Ideal  S  tructure  of  Ice  and 
Electrical  Properties  of  Ice  (Petrenko  1993a,b),  many 
important  properties  arise  from  the  disorder  of  the 
protons  in  the  structure.  This  allows  molecules  to 
rotate,  giving  rise,  for  example,  to  dielectric  polar¬ 
ization,  but  such  rotations  cannot  occur  with¬ 
out  the  aid  of  special  kinds  of  point  defects  in 
the  protonic  subsystem.  The  structure  of  ice  is 
shown  in  Figure  4a.  There  are  two  possible  posi¬ 
tions  for  protons  on  each  bond,  but  in  the  perfect 
structure  the  occupation  of  these  sites  is  governed 
by  the  so-called  Bemal-Fowler  rules  or  simply  the 
"ice  rules": 


1.  Two  protons  are  located  near  each  oxygen 
atom. 

2.  There  is  one  proton  per  hydrogen  bond. 

The  first  rule  ensures  electroneutrality  of  the  water 
molecules,  and  the  second  ensures  the  perfection 
of  the  hydrogen  bonds.  Violations  of  these  rules  at 
particular  bonds  or  molecules  result  in  the  forma¬ 
tion  of  the  four  protonic  defects  shown  in  Figure 
4b.  These  are  the  H30+  and  OH-  ions  and  the 
Bjerrum  L-  and  D-defects,  which  are  bonds  with 
no  protons  or  two  protons  respectively  (Bjerrum 
1951).  An  ion  moves  from  one  site  to  another 
merely  by  the  transfer  of  a  proton  from  one  end  of 
a  bond  to  another,  and  a  Bjerrum  defect  moves  by 
the  rotation  of  a  water  molecule. 

The  structure  of  the  hydrogen  bonds,  along 
which  the  protons  move,  constitutes  a  three-di¬ 
mensional  grid.  The  transfer  along  this  grid  of  any 
of  the  above  defects  reorients  water  molecules  in 
such  a  fashion  that  a  second  defect  of  the  same  type 
cannot  pass  along  the  same  path  in  the  same 
direction. 

To  better  understand  what  happens  in  the  pro¬ 
cess  of  ionic  and  Bjerrum  defect  transfer,  let  us 
consider  the  track  of  an  arbitrary  charge  carrier,  a 
positive  H30+  ion,  for  instance.  When  we  are 
speaking  about  ionic  or  Bjerrum  defect  transfer,  it 
is  important  to  realize  that  oxygen  atoms  do  not 
"hop"  from  one  lattice  site  to  the  next  one.  Only 
protons  move,  with  a  subsequent  redistribution  of 
the  electronic  density  on  the  neighboring  mol¬ 
ecules.  Thus,  when  speaking  about  H30+  motion, 
we  have  in  mind  that  a  characteristic  distribution 
of  three  protons  near  an  oxygen  atom  shifts  in 
space. 

When  a  proton  jumps  along  a  hydrogen  bond 
from  one  water  molecule  to  another,  a  pair  of 
H30+  and  OH-  ions  is  produced,  as  shown  in 
Figure  4b.  Formation  of  this  ion  pair  and  further 
separation  of  the  ions  require  considerable  energy 
of  about  1  eV  (see  the  review  of  experimental 
results  in  Petrenko  1993a).  After  an  ion  is  formed, 
it  can  migrate  from  one  site  to  another  through 
proton  jumping  along  hydrogen  bonds.  This  pro- 
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a.  Defect  free  structure  of  ice  Ift. 


Figure  4.  Ice  structures. 


cess  is  illustrated  schematically  for  the  straight¬ 
ened  out  chain  of  molecules  shown  in  Figures  5a 
andb. 

When  a  proton  jumps  from  one  bond  to  an¬ 
other,  two  neighboring  defective  hydrogen  bonds 
are  produced.  These  are  the  so-called  D-  and  L- 
defects.  Since  a  proton  has  a  positive  electrical 
charge,  it  is  not  difficult  to  see  that  a  D-defect  is 
positively  charged  and  an  L-defect  is  negatively 
charged.  Formation  and  separation  of  a  pair  of 
Bjerrum  defects  require  energy  of  0.68  eV.  Already 
formed  L-  and  D-defects  can  migrate  in  the  ice 
lattice  by  protons  jumping  between  neighboring 


hydrogen  bonds  as  illustrated  for  a  D-defect  in 
Figures  5c  and  d. 

An  H30+  ion  migrates  from  left  to  right  along  a 
water  molecule  chain  by  consecutive  jumps  of 
protons  along  hydrogen  bonds.  Therefore,  if  in  the 
beginning  the  protons  on  all  bonds  in  the  chain 
were  located  in  the  left-hand  positions  (Fig.  5a), 
after  an  H3O  ion  has  passed,  all  protons  will  be  in 
the  right-hand  positions.  This  path  is  now  closed 
for  proton  transfer  in  the  same  direction.  The 
reader  can  easily  make  sure  that  the  migration  of 
an  OH-  ion  from  right  to  the  left  would  result  in 
exactly  the  same  rearrangement  of  this  molecular 


a.  Initial  arrangement  ofH$0+  ion  and  water  molecules  in  the  chain  before  the  ion  moves 
from  the  left  to  the  right.  The  arrows  indicate  successive  proton  hops  (along  the  bonds)  by 
means  of  which  the  ion  moves. 


b.  Chain  after  the  ion  has  passed.  Obviously,  no  further  H$0+  ion  can  go  along  the  same 
chain  in  the  same  direction. 


c.  D-defect  motion  from  left  to  right  along  the  chain  as  in  Figure  5b.  The  arrows  indicate 
successive  hops  of  protons  from  one  bond  to  next  one. 


d.  Chain  after  the  D-defect  has  passed.  Comparison  with  Figure  5a  shows  that  motion  ofD- 
defect  restored  positions  of  the  protons  along  the  chain.  It  is  ready  (or  "open")  now  for  the 
passage  of  another  H3O+  ion. 

Figure  5.  Schematic  representation  ofHzO*  ion  (a,  b)  and  D-defect  (c,  d)  motion  along  a 
chain  of  hydrogen  bonds  in  ice. 


chain.  Since  H3O  and  OH-  ions  drift  in  opposite 
directions  in  an  external  electric  field,  we  come  to 
the  conclusion  that  both  types  of  ions  drifting  in  a 
given  field  "block"  molecular  tracks  in  the  same 
manner. 

If  we  now  consider  the  migration  in  the  same 
field  of  a  positively  charged  D-defect  along  the 
same  chain  (Fig.  5c)  "blocked"  by  an  H30+  ion 
passing  along  it,  it  turns  out  that  the  D-defect 
"unblocks"  the  chain,  restoring  the  initial  configu¬ 
ration  to  that  of  Figure  5d,  which  is  the  same  as 
Figure  5a.  In  the  same  fashion  the  migration  of  a 


negatively  charged  L-defect  "unblocks"  the  chain 
if  it  migrates  from  right  to  left. 

Ions  and  Bjerrum  defects  migrating  through  ice 
reorient  water  molecules,  altering  conditions  for 
their  transfer.  One  can  say  that  in  ice  protonic 
defects  interact  in  a  very  special  way:  via  those 
changes  in  the  molecule  configuration  that  they 
produce  during  their  transfer.  In  addition  the  more 
usual  interactions  between  defects  also  take  place. 
First,  protonic  defects  are  electrically  charged, 
which  results  in  a  strong  Coulomb  interaction,  and 
second,  as  we  shall  see  later  (in  the  Activation 
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Volume  of  Protonic  Defects  section),  the  presence  of 
a  defect  introduces  a  lattice  deformation  that  can 
lead  to  interactions  via  elastic  stresses. 

In  thermal  equilibrium,  which  is  characterized 
by  maximum  entropy  S,  water  molecules  are  ori¬ 
ented  randomly,  subject  to  the  ice  rules.  Any  flux 
of  defects  of  the  same  type  (HiO*  ions,  for  ex¬ 
ample)  j  j  will  increase  the  number  of  molecules 
oriented  along  the  direction  of  the  flux.  The  new 
configuration  of  water  molecules  is  more  ordered 
and  therefore  has  smaller  entropy.  These  ideas 
were  formulated  mathematically  by  Jaccard  (1959, 
1964)  and  independently  by  Onsager  and  Depuis 
(I960).  They  introduced  the  configurational  vector 
£2  describing  changes  occurring  in  ice  because  of 
protonic  defect  transfer.  If  at  time  t  =  0  ice  is 
unpolarized  and  in  thermal  equilibrium,  then 

^  rt 

Cl  =  {ji  -  ]2  -;-3  +/4  )  dt  (23) 

.  0 

where  ji  ( i  =  1,2,3, 4)  denote  the  fluxes  of  H3O, 
OH-  and  D-  and  L-defects  respectively.  This  nota¬ 
tion  of  protonic  defects  is  generally  accepted  and 
we  shall  use  it  hereafter.  If  we  introduce  the  quan¬ 
tities  r|,-  given  by 

r\i  =  1,-1 -1,1  for  i  =  1 ,2,3,4,  (24) 

then  eq  23  can  be  simplified  to 


The  choice  of  +  or  -  before  ji  in  eq  23  accounts 
for  different  effects  produced  by  defects  on  the 
reorientation  of  water  molecules.  Thus,  for  ex¬ 
ample,  as  we  have  seen  in  Figure  5,  the  effect  of  an 
H3O+  ion  ( i  =  1)  is  opposite  to  the  effect  of  a  D- 
defect  ( i  -  3). 

Configurational  entropy  Sc  must  somehow  be 
related  to  £2  .  For  small  decrements  of  entropy  we 
can  use  the  first  term  of  a  Taylor  series 

A  Scoc-Q2.  (26) 

There  is  no  linear  term  because  ASC  cannot 
depend  on  the  sign  of  the  fluxes  producing  it.  For 
the  time  derivative  of  the  entropy  we  have 

Soc-Q  Q  . 


The  proportionality  factor  in  eq  27,  which  we 
denote  as  <h/T,  can  in  principle  be  both  calculated 
theoretically  and  determined  experimentally.  At¬ 
tempts  to  calculate  O  theoretically  (Jaccard  1959, 
1964;  Hubmann  1979a,b)  yield  only  the  right  order 
of  magnitude  of  this  quantity,  and  the  most  reli¬ 
able  value  is  that  found  in  an  analysis  of  experi¬ 
mental  data  by  Hubmann  (1979a) 

j  =  3.85  kBr m  (28) 

where  rDO  =  2.76  A  is  the  oxygen-oxygen  distance 
in  ice. 

According  to  Ryzhkin*  the  difference  between 
the  theoretical  and  experimental  values  of  O/T  is 
caused  by  taking  into  account  reorientation  of 
molecules,  not  only  in  this  water  molecule  chain, 
but  also  reorientation  of  neighboring  molecules 
when  considering  a  defect  transfer  along  a  chain. 
However,  such  calculations  are  quite  cumbersome 
and  have  not  been  made  yet.  Accordingly,  we  shall 
use  the  value  of  O  found  experimentally. 

Charge  carriers  in  ordinary  ice  Ih  are,  under 
usual  conditions,  one  of  these  four  types  of  pro¬ 
tonic  defects.  A  detailed  discussion  on  the  nature 
of  electric  charges  in  ice  can  be  found  in  the  report 
by  Petrenko  (1993a).  "Usual  conditions"  are  un¬ 
derstood  here  as  an  absence  of  excitation  by  direct 
ultraviolet  light,  electron  beams,  y  rays,  etc. 

The  electric  current  density  J  in  ice  can  now  be 
written 

7=  <»> 

1=1 

where  et  are  the  effective  charges  carried  by  the 
defects.  Since  H3O  and  OH~  ions  and  L-  and  D- 
defects  are  created  and  annihilated  in  pairs,  it  is 
obvious  that 

ei  =  —62  and  e$  =  -e\  .  (30) 

If  we  place  an  extra  "foreign"  proton  into  an  oth¬ 
erwise  normal  ice  lattice,  we  will  inevitably  pro¬ 
duce  two  carriers  at  once:  an  H30+  ion  and  a  D- 
defect,  as  shown  in  Figure  6.  Hence,  the  sum  of 
these  charges  equals  a  proton  charge  e 

e\  +  es  =  e .  (31) 

‘Personal  communication  with  Dr.  Ivan  Ryzhkin  during  the 
International  Symposium  on  Physics  and  Chemistry  of  Ice, 
September  1991,  Sapporo,  Japan. 


(27) 
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Figure  6.  Creation  of  two  point  defects  (H^O+  ion  and  D-defect)  when  a  proton  is  placed 
into  defect-free  ice  lattice. 


To  find  out  how  the  elementary  protonic  charge  is 
shared  between  the  ion  and  the  D-defect,  we  must 
consider  both  the  fraction  of  the  intermolecular 
distance  that  the  proton  travels  in  the  case  of  ionic 
or  Bjerrum  defect  transfer  and  the  associated  dis¬ 
placement  of  electronic  charge  density  when  the 
defect  moves  along  the  chain. 

The  most  strict  theoretical  analysis  of  this  prob¬ 
lem  was  conducted  by  Scheiner  and  Nagle  (1983) 
who  found 

e3  =  (0.36  ±  0.03)e .  (32) 

The  other  way  to  find  e\  and  e3  is  by  experiment. 
In  ice  the  charge  of  the  predominant  defect  is 
incorporated  in  the  formula  for  static  dielectric 
permittivity  £$  (see  Petrenko  1993a).  Using  this  fact 
and  analysis  of  experimental  data  of  measure¬ 
ments  of  £s  for  specially  doped  ices  (in  such  a 
fashion  that  either  ions  or  Bjerrum  defects  pre¬ 
dominated),  Hubmann  (1979a,b)  found 

e3  =  (0.38  ±  0.01  )e  (33) 

which  is  remarkably  close  to  the  calculations  of 
Scheiner  and  Nagle.  Other  recent  experimental 
results  for  e\  and  e3  are  given  in  Table  4.  These 
results  are  close  to  those  represented  in  eq  32 
and  33. 


The  rate  of  heat  generation  in  a  unit  volume 
Q  by  current  flowing  through  ice  is  determined 
by  two  effects.  The  first  is  simply  the  Joule  heat 
E  J ,  and  the  second  is  the  change  in  the  configura¬ 
tional  entropy  TSC  Accordingly 


Q  =  e7-4>Hd. 

(34) 

Using  eq  25  and  29  we  obtain 

Q=  I  M-nfOa)  h. 

(35) 

i  =  1 


As  can  be  seen  from  eq  35,  the  effective  force  acting 
upon  charge  carriers  in  ice  consists  of  two  terms,  as 
opposed  to  ordinary  conductors 


Table  4.  Ratios  of  effective  electric  charge  of 
H30+  ion  (ei)  and  of  D-defect  (e3)  to  proton 
charge  (e). 

ei/e  e3/e  Reference 

0.61  ±  0.02  0.39  ±  0.01  Zaretskii  et  al.  (1988) 

0.62  ±  0.01  0.38  +  0.01  Hubmann  (1979a) 

0.626  0.374  Takei  and  Maeno  (1987) 

0.64  ±  0.03  0.36  ±  0.03  Scheiner  and  Nagle  (1983) 
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(41) 


Fj=  e,E- rj,4>n  . 


(36)  «D  =  «L- 


Accordingly,  for  the  charge  carrier  fluxes  we  have 
j  i  =  iei  E-  tl,4>  D )  n,  g,  le,  (37) 


Using  a  method  similar  to  the  one  seen  above, 
Ryzhkin  found  that  in  the  presence  of  both  types  of 
defects 


where  g;  is  the  mobility  (i.e.,  velocity  per  unit 
electric  field)  of  defect  i.  If  the  distribution  of 
protonic  defects  in  ice  is  not  uniform,  we  have  to 
add  the  diffusion  flux  to  eq  37,  giving 

j,  =  (c,£-t|,4>q)  /i,-g,7e,— Dj-gradn,-  (38) 

in  which  the  diffusion  coefficient  D,  is  equal  to 
kuTg/le,- 1.  The  four  equations  (eq  38,  i  =  1, 2,  3, 4) 
are  the  most  general  equations  describing  charge 
carrier  transfer  in  ice  and  can  account  for  all  elec¬ 
trical  properties  of  ice  known  today  (Petrenko 
1993a). 


Statistics  of  protonic  point  defects  in  ice 

To  calculate  the  equilibrium  concentration  of 
protonic  point  defects,  we  can  follow  the  same 
general  procedure  that  was  described  in  the  Gen¬ 
eral  Consideration  of  Point  Defects  section.  For  ice 
with  ionic  defects,  such  a  calculation  was  made  by 
Ryzhkin  (1985),  who  carefully  calculated  the  en¬ 
tropy  of  ice  with  H3O  and  OH-  ions.  Minimizing 
the  free  energy  F,  Ryzhkin  then  obtained  for  the 
equilibrium  concentration  of  ions 


— (5 —  =  2.  exp  ( — 
N-2n,  3  \  2k qT 


(39) 


where  Eaj  is  the  formation  energy  of  a  pair  of  ions, 
and  since  Ea;  »  2ku,T 


As  we  can  see,  the  specific  structure  of  ice  mani¬ 
fested  itself  only  in  the  appearance  in  eq  40  of  a 
factor  of  2/3  before  the  exponent.  We  can  interpret 
this  factor  in  simple  terms  as  the  ratio  of  the 
number  of  possible  orientations  of  an  ion  (four 
available)  to  the  number  of  possible  orientations  of 
a  water  molecule  (which  is  six).  The  factor  1  /2  is  in 
the  exponent  because  ions  are  created  in  pairs. 

The  general  case  of  the  presence  of  ions  and 
Bjerrum  defects  was  also  considered  by  Ryzhkin 
(unpublished  data,  1985).  Let  us  denote  the  D- 
defect  concentration  as  and  the  concentration 
of  L-defects  as  n  1.  If  D-  and  L-defects  are  formed 
by  thermal  fluctuations,  i.e.,  in  pairs,  then  obvi¬ 
ously 


S/(*bN)  =  -2  *  ln(x)  -  (1-2  x)  ln(2  (1-2  x)/3) 
-  4  y  ln(2  y)  -  2  (1-2  y)  ln(l-2  y)  (42) 


where 

x  =  n  JN  and  y  =  nn/N  .  (43) 

Minimizing  the  free  energy  of  ice  with  the  entropy 
from  eq  42,  we  obtain  for  the  Bjerrum  defect  con¬ 
centration 


«D="L=«H2oexp[-^;J.  (44) 

Protonic  defects  can  also  be  introduced  into  ice 
by  doping  with  certain  impurities.  This  case  will 
be  considered  in  the  Impurities  section. 

Electrical  properties  of  protonic  defects 

As  we  have  indicated  previously,  protonic  de¬ 
fects  are  charge  carriers  in  ice  and  hence  determine 
its  electrical  properties.  A  detailed  description  of 
the  electrical  properties  of  ice  and  experimental 
techniques  for  their  investigation  can  be  found  in 
the  recent  review  by  Petrenko  (1993a).  Here,  we 
shall  only  point  out  what  information  concerning 
protonic  defects  can  be  obtained  from  analysis  of 
the  electrical  properties  of  ice. 

When  an  alternating  electric  field  with  angular 
frequency  co  is  applied  to  ice,  in  addition  to  the 
drift  current  (eq  29\  J  <*,  we  have  to  add  the  dis¬ 
placement  current  J  dis,  so  that  the  total  current  is 


J  ~~  J  dis  +  J  dr 

(45) 

-V  3F. 

I  dis  ~  £o£oo”^7 
of 

(46) 

where  £0  is  the  dielectric  permittivity  of  free  space 
and  £00  is  the  high-frequency  dielectric  permittiv¬ 
ity  of  ice  arising  from  its  electronic  polarizability; 
&o  =  3.2. 

A  complete  description  of  the  electrical  proper¬ 
ties  of  ice,  i.e.,  its  electrical  conductivity  o(c 0)  and 
dielectric  permittivity  £(co),  can  be  found  in  the 
system  of  eq  25, 29, 45  and  46  (see  Petrenko  1993a, 
for  instance).  In  the  absence  of  interfaces  and  con¬ 
centration  gradients,  the  answer  can  be  expressed 
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as  follows 


(55) 


cr(co)  =  <rs  + 


(cr=o-  Os)(co  xDf 
1  +  (to  to)2 


e  (to)  =  £m  + 


1  +  (<o  td)2 


(47) 

(48) 


Here,  we  have  introduced  new  terms 
4 

Coo  =  X  °i  with  °i  =  ni  M  hi 
!=1 

e2/(Ts  =  _£ !ii_  +  ■  ...g.3.2_ 

Cl  +  <72  O3  +  O4 

7D_1  =  C0D  =  O 


Cl  +  Og  +  c 3  +  04 


eiz 


£s  —  £<»  + 


Cl  +  O 2  C3  +  <74 


ei 


e3 


j  Oi  +  <72  +  <73  +  er4\2 


ei 


(49) 

(50) 

(51) 


(52) 


By  varying  the  temperature  or  the  type  and 
concentration  of  dopants,  or  both  (see  the  Impuri¬ 
ties  section),  any  one  of  the  four  0;  terms  can  be 
made  dominant.  Then  eq  52  is  simplified  to  yield 

£s  ®  £00  4 — ,  (53) 

£0$ 

Measurements  of  £5  for  doped  ice  in  which  ionic 
defects  and  Bjerrum  defects  dominate  the  electri¬ 
cal  properties  thus  permit  both  the  charges  e;  and 
<5  to  be  determined.  This  is  the  basis  of  the  results 
quoted  in  Table  4. 

As  we  can  see  from  eq  49  and  50,  <7„  is  deter¬ 
mined  by  the  largest  0,-,  and  0S  is  determined  by 
the  lower  of  (01  +  02)  and  (03  +  04)  .  In  pure 
monocrystalline  ice  (03  +  04) » (01  +  02)  and  04  »  03. 
Consequently,  the  measurement  of  0M  as  a  func¬ 
tion  of  temperature  allows  us  to  determine  the 
activation  energy  of  the  partial  conductivity  04  of 
L-defects,  ECT4 

Coo  ~  04  =  040  exp  l^ij .  (54) 

Since  the  conductivity  of  Bjerrum  defects  greatly 
exceeds  (usually  by  two  or  three  orders  of  magni¬ 
tude)  the  ionic  conductivity,  and  the  mobility  of 
H30+  ions  exceeds  the  mobility  of  OH-  ions,  we 
have 


Therefore,  from  the  temperature  dependence  os(T), 
it  is  possible  to  calculate  the  conductance  activa¬ 
tion  energy  E01  of  H3O  ions.  Since  conductivities 
are  the  products  of  concentrations  and  mobilities, 
the  activation  energies  Eai  obtained  are,  generally 
speaking,  the  sums  of  formation  energies  (Eaj/2) 
and  activation  energies  of  migration  of  these  de¬ 
fects  Emj.  To  separate  out  (Eaj/2)and  Emi-,  usually  a 
known  concentration  of  specific  defects  is  achieved 
by  some  technique  and  then  the  0,  are  measured. 
The  known  concentration  of  the  defects  can  be 
created  with  proton  injection  (Petrenko  et  al.  1983), 
a  field  effect  transistor  (Petrenko  and  Maeno  1987) 
or  with  doping  of  ice  (see  the  Impurities  section). 
Thus,  in  many  cases  it  is  possible  to  determine 
separately  concentrations  «,•  and  mobilities  p,  of 
protonic  defects. 

The  most  well-established  conductivity  param¬ 
eters  of  pure  monocrystalline  ice  are  the  high- 
frequency  conductivity  0 u  and  its  activation  en¬ 
ergy  Ecoo.  For  T  >  -50°C  these  depend  only  on  the 
L-defects.  These  quantities  can  be  well  reproduced 
because,  first,  high-frequency  measurements  do 
not  depend  upon  the  type  of  electrode  used  and, 
second,  because  of  the  comparative  ease  with  which 
specimens  can  be  obtained  that  are  pure  enough 
for  the  L-defect  concentration  to  be  intrinsic.  The 
same  holds  for  td  and  its  activation  energy  ET, 
since  according  to  eq  49-51,  td-1  ^  0*.  when  0*, » 
0S  ( i.e.,  when  the  majority  type  of  charge  carrier  is 
clearly  distinguished).  These  four  values  at  T  = 
-10°C  are  as  follows 

ECToo=  (0.58  ±  0.03)eV 

Et  =  (0.58  ±  0.01)eV 

0oo  =  (1.8  ±  0.2)10-5  O-l  m-1 

td  ~  5  x  10-5  s  .  (56) 

The  Debye  dispersion  shown  in  Figure  7  is  repro¬ 
duced  very  well  in  good  quality  single  crystals: 
pure,  homogeneous  and  annealed  (for  relaxation 
of  internal  stresses). 

Much  greater  scatter  is  observed  in  measure¬ 
ments  of  the  static  conductivity  0S  of  undoped  ice 
and  its  activation  energy  Effs.  Values  of  Effs  have 
been  found  to  lie  between  0  and  0.7  eV  and  0S 
(-10°C)  varies  from  10 T6  down  to  6  x  10~10  D_1 
m_1 .  This  variation  is  explained  by  the  difficulty  of 
purifying  ice  to  a  level  at  which  intrinsic  ions 
dominate  over  ions  supplied  by  impurities.  That  is 
why  von  Hippel  et  al.  (1973)  concluded  that  there 


13 


-4 


-5 


-6 


-7 


-8 


0 

o> 

o 


Figure  7.  Typical  frequency  dependence  of  permittivity  e  and  of  conductivity  a  of 
pure  single  crystalline  ice  at  T  =  -10C°  (after  Petrenko  1993a). 


are  no  intrinsic  ions  in  ice  and  that  all  ions  origi¬ 
nate  from  impurities.  If  we  refer  to  Table  I  of  their 
paper  (von  Hippel  et  al.  1971)  for  the  value  of  crs 
(2.5  x  10~8  £2_1  m_1)  and  compare  it  with  the 
minimum  observed  value  of  6.4  x  10-10  Q_1  nr1, 
then  we  might  agree  with  them  that  their  ions  were 
of  impurity  origin.  The  question,  however,  re¬ 
mains  open:  were  the  specimens  in  the  experi¬ 
ments  of  Worz  and  Cole  (1969)  or  Petrenko  et  al. 
(1983)  pure  enough  that  their  conductivity  may  be 
considered  intrinsic?  We  observed  the  following 
general  trend:  the  purer  the  ice  is,  the  smaller  is  cts 
and  the  greater  is  Eas.  Therefore,  the  value  of  the 
activation  energy  of  creation  of  an  ion  pair  Eai  = 
0.96  eV  (Eigen  and  De  Maeyer  1956,  1957,  1958; 
Eigen  et  al.  1964),  for  a  long  time  accepted  as  the 
true  value,  should  be  treated  as  the  lowest  possible 
limit  of  this  quantity.  As  the  current  upper  limit  we 
have  Ea,  =  2Eas  =  1.4  eV  from  the  work  of  Petrenko 
et  al.  (1983).  The  above  discussion  implies  that  ion 
mobility  almost  does  not  depend  on  temperature, 
that  is,  the  ion  motion  is  not  a  thermoactivated 
process.  That  suggestion  was  proven  by  experi¬ 
ment  (see,  e.g.,  Kunst  and  Warman  1983). 

The  determination  of  charge  carrier  mobilities 
in  ice  is  an  extremely  difficult  task,  since,  because 
of  small  values  of  p,-  and  the  "jump"  mechanism  of 
motion,  the  regular  technique  for  measuring  the 
Hall  effect  is  not  applicable  in  this  case.  (The  Hall 
effect  is  discussed  in  the  future  report  on  the 


Surface  of  Ice.)  Hopes  of  using  "saturation  cur¬ 
rents"  for  determining  p,  are  not  justified  because 
of  reasons  pointed  out  by  Maidique  et  al.  (1971)  in 
their  critique  of  Eigen's  work,  and  also  because 
when  high  voltage  is  applied  to  a  system  with  ice/ 
metal  interfaces,  various  physical  phenomena  are 
observed  that  make  the  interpretation  of  results 
more  complicated  (see  Part  II  of  Petrenko  1993a). 

Nonetheless,  several  new  techniques  for  deter¬ 
mining  charge  carrier  mobilities  in  ice  have  been 
developed  and  employed  during  the  last  25  years. 
These  are:  the  technique  that  measures  stationary 
and  transient  injection  currents  (Eckener  et  al. 
1973,  Petrenko  et  al.  1983,  Petrenko  and  Ryzhkin 
1984b);  the  ice  field  effect  transistor  (Petrenko  and 
Maeno  1987);  and  the  technique  that  measures 
low-frequency  dispersion  (Petrenko  and  Ryzhkin 
1984a,  Zaretskii  et  al.  1988,  Zaretskii  1991).  The 
idea  of  determining  p,  and  «,  from  the  parameters 
of  low-frequency  dispersions  was  presented  by 
Petrenko  (1993a).  Other  techniques  mentioned 
above  are  described  and  discussed  in  Part  II  of  the 
same  review.  The  results  of  determining  charge 
carrier  mobilities  in  ice  are  summarized  in  Table  5. 
At  first  glance  only  the  data  on  L-defect  mobilities 
seem  to  be  satisfactory.  At  T  =  -10°C,  P4  =  (2  -  5)  x 
10-8  m2  V-1  s-1.  This  mobility  drops  exponentially 
with  decreasing  temperature,  the  activation  en¬ 
ergy  being  0.19  -  0.235  eV  (see  Table  6).  There  is 
still  no  reliable  information  about  D-defect  mobil- 
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Table  5.  Mobilities  of  charge  carriers  in  ice. 


T 

(°C) 

H3O+ 

P-i 

(m^/V  s) 

OH- 

(nfi/Vs) 

D 

in 

(mWs) 

L 

\U 

{mW  s) 

Method  used 

Reference 

-13...  -36 

(1.1  ±0.1)x  10-7 

— 

— 

— 

Analysis  of 
conductivity  decay 

Kunst  and  Warman  (1983) 

-10 

2.5  x  10-8 

— 

— 

— 

Os  and  Zias 

measurements 

Worz  and  Cole  (1969) 

O 

*? 

O 

5  x  10-7 

— 

— 

— 

Saturation  currents 

Maidique  et  al.  (1970) 

-123 

_ 

<4xl(H 

— 

— 

Doping  with  KOH 

Zaretskii  et  al.  (1988) 

-5...  -40 

10-9-  10-7 

— 

— 

— 

Proton  injection 

Petrenko  et  al.  (1983) 

0 

2.7  x  10-8 

— 

— 

5  x  10-8 

Doping  with  HF 

Camplin  et  al.  (1978) 

-10 

_ 

— 

— 

1.8x10-8 

Doping  with  HF 

Jaccard  (1959) 

— 

— 

>2x10-10 

— 

— 

Doping  with  KOH 

Howe  and  Whitworth 
(1989) 

-10 

(2.4  ±  1.6)  x  10-7 

— 

— 

— 

Saturation  currents 

Bullemer  et  al.  (1969) 

-10 

3.5  x  10-8 

— 

— 

— 

HF  doping 

Camplin  and  Glen  (1973) 

-33.1 

9.2  x  10-8 

2.7  x  10-8 

— 

— 

Ice  field  effect  transistor 

Petrenko  and  Maeno 
(1987) 

-145...  -178 

(9  ±  1)  x  10-8 

— 

— 

— 

Transient  injection  of 
protons 

Eckener  et  al.  (1973) 

-10 

=  10-6 

— 

— 

1.7  x  10-8 

Low-frequency  dispersions 

Zaretskii  et  al.  (1987a) 

-33 

(6.7  ±  0.8)  x  10-7 

— 

— 

— 

Low-frequency  dispersions 

Zaretskii  (1991) 

-20 

6  x  10-8 

— 

— 

— 

Proton  injection 

Petrenko  and  Ryzhkin 
(1984c) 

-10 

— 

2.4  x  10-8 

— 

— 

Recombination  injection 

Petrenko  and  Chesnakov 
(1990c) 

-10 

(9  ±  3)  x  10-8 

7  x  10-9 

— 

_ 

Analysis  of  data 

Chen  et  al.  (1974) 

ity.  There  are  only  indirect  evidences  that  14.3  <  (14. 
For  example,  Devlin  (1992),  citing  spectroscopic 
data,  indicated  that,  in  the  region  at  150  K,  p.3  is 
orders  of  magnitude  less  than  (J.4- 

Data  on  ionic  mobilities  are  distinguished  at 
first  glance  by  their  large  scatter  (from  5  x  10-7 
down  to  2.5  x  10-8  for  HaO  and  from  4  x  10-7 
down  to  2  x  HH°  for  OH-) .  However,  if  we  discard 
the  large  values  obtained  by  the  now-discredited 
method  of  saturation  currents  and  the  values  ob¬ 
tained  from  measurements  of  low-frequency  dis¬ 
persions  (which  still  have  low  accuracy),  it  ap¬ 
pears  that  the  positive  ion  mobility  lies  within  the 
limits  of  2.7  x  lO8  and  1.1  x  10~7  m2/V  s — not  bad 
at  all!  The  authors  of  this  review  prefer  the  data 
obtained  in  transient  injection  experiments 
(Eckener  et  al.  1973),  in  which  a  minimum  of 
assumptions  was  used.  The  same  HsO+  ion  mobil¬ 
ity  was  found  by  Petrenko  andMaeno  (1987)  using 
field  effect  transistors.  So  it  is  likely  that  pi  =  10-7 
m2  V-1  s-1.  It  seems  well  established  that  both 
H3O+  and  OH-  ions  move  in  ice  with  zero  activa¬ 
tion  energy. 

Concluding  this  section  about  point  defects' 
mobilities,  we  can  now  explain  the  at  first  slightly 


surprising  fact  that  the  self-diffusion  coefficients 
of  oxygen  and  hydrogen  in  ice  are  equal  (Table  2). 
This  is  unexpected  because  the  oxygen  can  diffuse 
only  as  a  part  of  a  water  molecule,  whereas  the 
hydrogens  can  also  move  by  ion  and  Bjerrum 
defect  migration.  Let  us  compare  the  absolute 
values  of  diffusion  coefficients  of  hydrogen  in 
these  two  processes.  If  the  hydrogens  diffuse  in  ice 
by  protonic  defect  transfer,  then  the  diffusion  co¬ 
efficient  Dh  is 

Dh  =  nipi  +  K2°2  (57) 

nH20 

where  n\i2  and  D\/2  denote  the  concentrations  and 
diffusion  coefficients  of  those  protonic  defects  that 
limit  the  transfer  of  protons  over  macroscopic 
distances.  In  pure  ice  these  defects  are  ions.  There¬ 
fore,  we  can  rewrite  eq  57  in  the  form 

Dh  =  ei(niPi  +  n2[i2)kBT  =  crs •  khT  (5g) 
^HzO  el  *  nH2o 

where  as  is  the  static  electric  conductivity  of  ice. 
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For  Dh  to  exceed  the  self-diffusion  coefficient  Ds, 
the  following  inequality  must  hold 

<7S  >  pse2”H2°  *  3.9  x  10-*  (Q  m)-1  .  (59) 

k#T 

In  eq  59  we  have  used  Ds  =  2  x  1(H5  m2  s_1  (T  = 
-10°C)  from  Table  2.  Since  in  pure  ice  at  -10°C  crs 
<  10~8  (ft  m)”1  (Petrenko  1993a),  the  diffusion  rate 
of  hydrogen  in  ice  is  faster  when  it  diffuses  as  a 
part  of  a  water  molecule  in  the  form  of  an  intersti¬ 
tial.  That  explains  the  equality  of  the  diffusion 
coefficients  of  hydrogen  and  oxygen.  However,  in 
ice  heavily  doped  with  HF  or  KOH,  <rs  would  be 
sufficiently  large  for  hydrogen  transfer  by  pro¬ 
tonic  defects  to  be  faster  than  that  by  interstitial 
diffusion  of  H2O  molecules.  There  are  no  experi¬ 
mental  data  on  self-diffusion  coefficients  in  doped 
ice. 

Activation  volume  of  protonic  defects 

The  formation  of  any  protonic  point  defect  re¬ 
sults  in  deformation  of  the  surrounding  lattice.  In 


the  case  of  Bjerrum  defects,  the  defective  hydro¬ 
gen  bonds  are  weaker  than  the  normal  ones,  and 
therefore  they  will  be  longer.  This  is  obvious  in 
particular  for  a  D-defect,  where  there  exists  an 
extremely  strong  repulsion  between  two  protons 
on  one  hydrogen  bond.  Accordingly,  one  can  as¬ 
sume  that  Bjerrum  defects  "expand"  ice,  increas¬ 
ing  its  volume.  If  ice  is  under  pressure  P  then  the 
formation  energy  of  the  2th  defect  Ea;-  is  increased 
by  the  amount  of  work  against  this  pressure 

Eaj(P)  =  Eai(0)  +  Pyaj  (60) 

where  ya;  is  the  so-called  activation  volume.  For 
the  defects  that  increase  volume,  ya,  is  positive  and 
their  formation  energy  is  raised  proportionally  to 
the  applied  pressure.  This  results  in  an  exponen¬ 
tially  decreasing  concentration  of  such  defects  with 
increasing  pressure  (Table  7). 

It  is  a  much  more  difficult  problem  to  predict 
even  the  sign  of  ya,  for  ions  in  ice.  It  was  believed 
initially  that,  owing  to  Coloumb  attraction  of  the 
dipole  moments  of  neighboring  water  molecules. 


Table  6.  Data  on  high-frequency  (a»)  and  static  (os)  conductivities  of  ice. 


Type  of  icc 

T 

(°C) 

(Too 

(ftm)-i 

as 

(Qm)-i 

Es» 

(eV) 

Eos 

(eV) 

4  zone  refinings 

-10 

— 

4  x  10-* 

_ 

0.54 

single  crystals 

— 

— 

2.5  x  10-8 
(0...-20°C) 

— 

0-0.06 

single  crystals 

-A  ...  -56 

— 

6.4  x  10-10 
(-lOX) 

0.53  ±  0.01 

0.70  ±  0.07 

single  crystals 

-10 

— 

10-7 

0.37 

— 

single  crystals 

— 

— 

0.624 

— 

single  crystals 

-  10 

1.6x10-5 

1.1  X  10-7 

0.575 

0.61 

pure  and  HF-doped 

single  crystals 

0 

— 

2.5  x  10-s 

— 

— 

single  crystals 

-10 

— 

(1.1  ±05)  x  10-* 

0.76  ±  0.06 

_ 

single  crystals 

— 

— 

2.3  x  10-8 

— 

0.49 

single  crystals 

-10 

1.8  x  10-5 

— 

0.61 

_ 

single  crystals 

-10 

— 

— 

— 

— 

HC1  doping 

£-39 

— 

— 

0.492  ±  0.025 

0 

10s  mol/L 

HC1  doping 

-39 

— 

— 

0.259  ±  0.003 

0.267  ±  0.009 

KOH  doping 

0.1  mol/L 

0...  -150 

2x10-5-101 

up  to  10-1 

0 

(TS-  70°C) 

0.85 

(T  <  -  70°C) 

(in  solution) 

KOH  doping 

-23...  -203 

-2x10-* 

1(H 

0 

0.4 

0.032  mol/L 

(-23°C) 

(-23°C) 

(T  >  -  100°C) 

(T  >  -  100°C) 

(in  solution) 
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ions  will  contract  the  lattice,  i.e.,  ya;  <  0.  Experi¬ 
mental  data  shown  in  Table  7  refuted  this  opinion, 
and  the  activation  volumes  of  both  H3O  and  OH~ 
ions  turned  out  to  be  positive.  We  can  try  to 
explain  this  by  the  fact  that  dipole  moments  of  all 
neighboring  molecules  in  the  ice  structure  cannot 
be  oriented  towards  the  ion.  This  would  result  in 
a  violation  of  the  Bemal-Fowler  rules.  Also,  if  a 
dipole  moment  is  directed  away  from  the  ion,  it 
will  be  repelled  rather  than  attracted,  increasing 
the  volume  around  the  defect.  Thus,  such  simple 
considerations  are  not  sufficient  to  estimate  ya!-  of 
ions,  and  ab  initio  calculations  taking  into  account 
the  permitted  configurations  of  water  molecule 
orientations  in  ice  are  required. 

For  a  defect  to  jump  from  one  molecule  to  the 
neighboring  one  (or  from  one  bond  to  another 
one),  an  additional  change  in  volume  ymj  may  be 
required.  This  is  rather  obvious,  since  the  energy 
of  the  transitional  state  of  a  defect  is  of  course  not 
lower  than  the  energy  of  the  equilibrium  state. 
Therefore,  when  considering  a  defect  transfer  un¬ 


der  pressure,  we  have  to  add  an  additional  term  to 
the  activation  energy  of  their  motion  EmI- 

Fmi (P)  =  Emi(0)  +  Ym i  P  ■  (61) 

In  other  words,  the  amplitude  of  the  potential 
relief  shown  in  Figure  1  varies  with  pressure. 
During  ion  transfer,  which,  according  to  experi¬ 
mental  data,  occurs  through  a  non-activational 
hop  of  protons  along  a  bond,  the  crucial  factor  may 
be  an  exponential  increase  in  the  probability  of 
such  a  transition  attributable  to  a  decrease  in  the 
length  of  the  hydrogen  bond  owing  to  the  applied 
pressure.  This  would  appear  as  a  negative  activa¬ 
tion  volume  ym,-. 

In  the  course  of  studies  on  the  dependence  of 
partial  conductivities  of  defects  a,  on  pressure,  it 
was  found  that 

Oi  =  \e\  wn  - 

eXp  |[£ai(Q)  +  Emi(Q)]  +  ^[y a;  +  Ymi]|  ^£) 


Table6(cont'd). 


Efli 

(eV) 

E „B 
(eV) 

Activation 
energy  of 
ion  mobility 
(eV) 

Activation 
energy  of 

D-  and  L-defect 
mobility 
(eV) 

T D 

(s) 

Et 

(eV) 

Reference 

1.08 

_ 

_ 

_ 

_ 

_ 

Wore  and  Cole  (1969) 

1.1 

— 

— 

— 

5.18  x  10-5 

0.598 

Maidique  et  al.  (1970) 

(-10°C) 

1.4 

— 

— 

— 

— 

— 

Petrenko  et  al.  (1983) 

_ 

_ 

_ 

— 

— 

- - 

Bullemer  and  Riehl  (1966) 

— 

— 

— 

— 

5  x  10-5 

— 

Caraplin  et  al.  (1978) 

1.22 

0.68 

0 

0.235 

2  x  lO-’ 

0.575 

Jaccard  (1959) 

(0°C) 

' 

— 

— 

— 

— 

1.8x10-5 

0.59 

von  Hippel  et  al.  (1971) 

Maidique  et  al.  (1971) 

— 

— 

— 

— 

— 

— 

Bullemer  et  al.  (1969) 

0.98 

— 

— 

— 

Maeno  (1973) 

— 

— 

— 

— 

4.5  x  10-5 

0.64 

Camplin  and  Glen  (1973) 

— 

— 

— 

— 

10-4 

0.57 

Zaretskii  et  al.  (1989) 

— 

— 

0 

— 

— 

0.567  ±  0.00 

Takei  and  Maeno  (1984) 

90.293  ±  0.034 

— 

0.79 

0 

0.19 

— 

— 

Takei  and  Maeno  (1987) 

— 

— 

0 

— 

2x10-5 

0 

Zaretskii  et  al.  (1988) 

(T<-  70°C) 

_ 

_ 

_ 

_ 

_ 

_ 

Howe  and  Whitworth 

0 

(1989) 
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Table  7.  Activation  volumes  of  protonic  point  defects  in  ice.  The  data  are  given  in  the  units  of  Vmol 
=  3.3  x  10-23  cm3. 


Ice  type 

( concentrations 
in  mole/L) 

ya 

( HjO *  ion ) 

y«i 

(OH  ion) 

y«) 

(D-defecI) 

y ^ 

(L-  defect) 

y 'at. 2  +  Ym/,2 

Y aJ.4  ±  Y 'm3, 4 

Reference 

pure  polycrystals 

— 

— 

— 

— 

-(0.55  ±0.15) 

0.145  ±0.015 

Chen  et  al.  (1965) 

pure  single  crystals 

— 

— 

— 

— 

-(0.18  ±  0.01) 

0.13  ±0.01 

Taubenberger  et  al. 
(1973) 

=  5  x  10-fi  HF 

— 

— 

— 

— 

0.091  ±0.024 

-(0.16  ±0.024) 

Taubenberger  et  al. 
(1973) 

=  3x  10-5  HF 

— 

— 

— 

— 

0 

-(0.36  ±  0.05) 

Taubenberger  et  al. 
(1973) 

pure  single  crystals 

— 

— 

— 

— 

-(3  ±  0.9) 

=  0.15 
(~25°C) 

Hubmann  (1978) 

10-4  NHj 

— 

— 

— 

— 

=  -0.7 

=  0.27 
(~25°C) 

Hubmann  (1978) 

3  x  10-4  hf 

— 

— 

— 

— 

=  -0.55 

=  0.15 
(-25°C) 

Hubmann  (1978) 

single  crystals 

10-4  HC1 

3.3  ±0.12 

— 

— 

1.33  ±0.3 

— 

— 

Evtushenko  and 
Petrenko  (1991) 

single  crystals 

10-4  HF 

3.9  ±  0.36 

— 

— 

1.0  ±0.39 

— 

— 

Evtushenko  and 
Petrenko  (1991) 

single  crystals 

10-4  NHj 

— 

0.82  ±0.39 

0.7  ±0.3 

— 

— 

— 

Evtushenko  and 
Petrenko  (1991) 

We  can  determine  only  the  "total"  activation  vol¬ 
ume  [yal-  +  Ymil-  Such  measurements  were  origi¬ 
nally  performed  by  Chen  et  al.  (1965)  on  pure 
polycrystalline  ice.  They  found  that  the  static  con¬ 
ductivity  cts  increases  and  the  high-frequency  con¬ 
ductivity  (To,,  decreases  with  increasing  pressure. 
Similar  measurements  on  monocrystalline  ice  were 
performed  later  by  Taubenberger  et  al.  (1973)  and 
Hubmann  (1978).  Taubenberger  et  al.  used  pure 
and  HF  doped  ice  and  Hubmann  used  pure  and 
NH3  and  HF  doped  ice.  The  idea  of  using  the 
doped  ice  was  to  produce  a  known  initial  concen¬ 
tration  of  various  defects  to  separate  out  ya;-  and 
ym,\  However,  the  interpretation  of  the  results 
obtained  from  the  doped  ice  still  remains  compli¬ 
cated  and  ambiguous,  since  pressure  can  affect  the 
dissociation  energies  of  protonic  defects  from  im¬ 
purities.  Table  7  contains  values  of  (ya; + ymI)  found 
in  these  works.  Here,  we  assume  that  and 
Ter1  are  proportional  to  the  partial  conductivity  of 
the  majority  charge  carrier  and  that  cts  is  propor¬ 
tional  to  the  partial  conductivity  of  the  minority 
charge  carrier,  in  conformity  with  eq  49-51 .  In  the 
analysis  of  data  presented  in  Table  7,  we  have  to 
take  into  account  that  in  the  case  of  doped  ice  Eai  is 
related,  obviously,  to  the  energy  of  separation  of  a 
defect  from  an  impurity. 


Evtushenko  and  Petrenko  (1991)  used  a  quite 
different  approach  to  determine  activation  vol¬ 
umes  of  protonic  defects  in  ice.  Since  the  formation 
energy  of  defects  Ea!-  depends  on  the  applied  pres¬ 
sure  P  according  to  eq  57,  in  the  presence  of  a 
pressure  gradient  there  is  a  force  F,  acting  upon  the 
defects 

=  -gra d(Ea, )  =  -ya;  grad(P) .  (63) 

In  their  work,  Evtushenko  and  Petrenko  intro¬ 
duced  a  pressure  gradient  by  bending  thin  ice  sam¬ 
ples  and  then  measured  the  electric  field  strength 
E  compensating  for  the  effect  of  the  force  F\ 

ei  E  =  -yai  grad(P) .  (64) 

By  varying  temperature  and  doping,  they  deter¬ 
mined  activation  volumes  of  all  four  types  of 
protonic  defects  in  ice  (seeTable  7).  The  advantage 
of  this  method  is  that  it  enables  us  to  separate  out 
Ea,  (i.e.,  the  additional  term  to  the  specific  volume 
of  a  defect)  from  the  sum  EaI  +  Em/.  The  results  pre¬ 
sented  in  the  table  are  expressed  in  terms  of  the 
volume  occupied  by  a  water  molecule  in  ice:  3.3  x 
IO23  cm3.  According  to  the  results  of  Evtushenko 
and  Petrenko,  all  protonic  defects  (except  H30+ 
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ion)  have  approximately  doubled  volume  com¬ 
pared  to  the  molecular  volume. 

Spectroscopic  study  of  protonic  defect 
migration — D-H  exchange 

An  original  and  informative  optical  method  of 
investigating  protonic  defects  in  ice  was  devel¬ 
oped  by  Devlin  (Ritzhaupt  and  Devlin  1977; 
Ritzhaupt  et  al.  1978,  1980;  Devlin  1990).  The 
method  is  based  on  the  optical  detection  of  hydro¬ 
gen  isotope  exchange  (H  and  D)  between  water 
molecules  in  ice.  Since  such  an  exchange  can  occur 
only  through  protonic  defect  transfer,  the  rate  of 
isotope  exchange  enables  us  to  estimate  the  diffu¬ 
sion  of  ions  and  Bjerrum  defects  in  ice.  Figure  8 


illustrates  the  essence  of  this  method.  If  we  place  a 
D2O  molecule  in  H2O  ice,  then,  because  of  the  trans¬ 
fer  of  ions  displacing  protons  and  deuterons  along 
hydrogen  bonds,  pairs  of  HOD  molecules  will  be 
formed  on  adjacent  sites.  We  shall  denote  these 
pairs  as  (HOD)2-  Hence,  the  rate  of  the  reaction 

D20  +  H20  ->  (HOD)2  (65) 

will  be  determined  by  the  number  of  intermolecu- 
lar  jumps  by  ions,  which  is  proportional  to  the 
ionic  conductivity  0iOn- 

For  two  deuterons  to  move  a  large  distance 
apart,  it  is  necessary  to  transfer  a  deuteron  from 
one  bond  to  another.  This  can  take  place  only 


Example  of  H,(t  ion 
and  L-  defect 


1  L-defect  migration  which 

H  separates  the  D  by  another 


oxygen. 


Figure  8.  Point  defects  in  an  ice  lattice,  showing  the  proton  hop  step  and  the  L-defect  turn  step 
of  proton  transport  and  isotopic  exchange. 
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through  Bjerrum  defect  migration.  Therefore,  the 
rate  of  the  reaction 

(HOD)2  HOD  +  HOD  (66) 

will  be  determined  by  the  number  of  elementary 
jumps  of  Bjerrum  defects  per  second,  which  is 
proportional  to  pb  «b  oc  ob- 

The  reactions  (eq  65-66)  can  be  observed  be¬ 
cause  the  infrared  absorption  spectra  of  D20, 
(HOD)2  and  isolated  HOD  are  different,  as  shown 
in  Figure  9.  Devlin  and  co-authors  prepared  H20 
ice  samples  with  a  small  addition  of  D20  by  vapor 
deposition  in  a  vacuum  at  125  K.  At  such  a  low 
temperature,  the  protonic  defect  transfer  is  negli¬ 
gibly  small  and  the  ice  obtained  has  the  cubic 
structure  Ic.  Varying  the  temperature  in  5-K  steps 
and  observing  the  kinetics  of  the  changes  in  the  IR 
spectra,  Devlin  and  co-authors  (Collier  et  al.  1984) 
found  that  in  their  samples  in  the  temperature 
range  from  130  to  150  K 


Figure  9.  Infrared  spectra  for  isotopomers  o/H20 
decoupled  in  H20  ice  at  90  K:  (a)  D20;  (b)  (HOD)2; 
(c)  HOD  (after  Devlin  1992). 


Eoi  =  0.41  eV 

and  (67) 

EoB  =  0.52  eV 

which,  although  close  to  the  values  determined 
from  electrical  conductivity  (see  Table  6),  are  still 
considerably  lower  than  the  conventional  values 
at  higher  temperatures.  This  discrepancy  may  be 
attributed  to  a  number  of  factors,  including  extrin¬ 
sic  defect  activity,  the  cubic  nature  of  the  ice  used, 
or  the  predicted  behavior  ignoring  shallow  proton 
trapping,*  or  all  three. 

Probably  the  major  significance  of  the  work  of 
Devlin  and  colleagues  is  that  they  confirmed  at  the 


*The  shallow  proton  traps  were  identified  with  the  nega¬ 
tively  charged  L-defects  (Kunst  and  Warman  1983).  An  asso¬ 
ciation  energy  of  a  proton  on  such  a  trap  was  found  as  =0.43  eV 
by  Woodbridge  and  Devlin  (1988)  and  =  0.19  eV  by  Kunst  and 
Warman  (1983). 


Figure  10.  Oxygen  atom  positions  in  ice  Ip. 
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molecular  level  the  two-step  character  of  proton 
transfer  in  ice,  and  consequently  the  necessity  for 
such  a  transfer  to  be  performed  by  two  groups  of 
very  different  defects:  ions  and  Bjerrum  defects.  In 
other  words,  the  work  of  this  group  provides  very 
strong  support  for  the  existence  of  both  types  of 
protonic  defects  in  ice.  Readers  interested  in  this 
question  can  find  details  in  the  review  by  Devlin 
(1990).  We  shall  also  return  to  this  topic  in  the 
report  Optical  Properties  of  Ice. 

Atomic  structure  of  protonic  defects 

The  atomic  structure  of  protonic  defects  de¬ 
picted  in  Figure  4b  is  purely  schematic.  The  intro¬ 
duction  of  a  defect  changes  considerably  the  forces 
acting  between  adjacent  atoms,  and  this  leads  to 
lattice  "relaxation"  around  the  defect.  One  piece  of 
evidence  for  such  relaxation  is  the  increased  spe¬ 
cific  volume  of  defects  as  compared  to  the  molecu¬ 
lar  volume  (see  the  Activation  Volume  of  Protonic 
Defects  section).  The  necessity  for  atomic  rear¬ 
rangement  of  the  ice  lattice  around  protonic  de¬ 
fects  can  be  demonstrated  by  estimation  of  the 
electrostatic  energy  of  interaction  between  two 
protons  on  the  same  hydrogen  bond  in  a  D-defect. 
If  we  do  not  take  relaxation  and  lattice  rearrange¬ 
ment  into  account,  then  the  two  protons  will  be  at 
a  distance  r  of  about  0.7  A  from  each  other.  Since  at 
such  small  distances  we  have  to  use  e  =  1,  the 
electrostatic  energy  of  such  a  D-defect  would  be 

AE  -  sL.  =  19  eV.  (68) 

Akteq, 

This  is  extremely  large  compared  to  the  observed 
formation  energy  of  0.68  eV  of  a  pair  of  L-  and  D- 
defects. 

Currently,  there  are  two  models  in  which  the 
energy  of  protonic  defects  is  lowered  by  rear¬ 
rangement  of  atoms  around  a  defect.  In  the  first 
model,  after  a  protonic  defect  is  formed  in  the 
simplified  fashion  shown  in  Figure  4b,  atoms  and 
molecules  are  allowed  to  shift  slightly  from  their 
initial  positions.  The  new  equilibrium  stable  con¬ 
figuration  is  the  one  with  the  lowest  energy.  This 
approach  to  the  calculation  of  energy  and  struc¬ 
ture  of  Bjerrum  defects  in  ice  was  used  by  Plummer 
(1987, 1992).  She  used  ab  initio  computer  simula¬ 
tions,  requiring  very  large  amounts  of  computer 
time,  and  because  of  this  only  clusters  comprising 
eight  water  molecules  with  a  Bjerrum  defect  in  the 
middle  could  be  simulated.  Two  ways  of  choosing 
such  clusters  with  "eclipsed"  and  "staggered" 
structure  are  shown  in  Figure  10.  The  new  inter¬ 
atomic  distances  and  energies  AE  for  the  formation 


Table  8.  Calculated  parameters  of  Bjerrum 
defects  (after  Plummer  1992). 


Parameter 

Eclipsed 

Staggered 

0-0  distance  (A) 

D-defect 

3.30 

2.72 

H-H  distance  (A) 

1.78 

2.23 

Creation  energy  (eV) 

0.417 

0.278 

0-0  distance  (A) 

L-defect 

3.18 

3.16 

Creation  energy  (eV) 

0.317 

0.343 

of  L-  and  D-defects  in  such  clusters  are  summa¬ 
rized  in  Table  8. 

Since  in  a  real  ice  lattice  Bjerrum  defects  can  be 
located  at  hydrogen  bonds  both  parallel  (eclipsed 
structure)  and  inclined  (staggered  structure)  to  the 
hexagonal  c-axis,  we  should  take  some  average 
between  the  eclipsed  and  staggered  values  for 
comparison  with  experimental  values,  taking  into 
account  that  there  are  three  times  more  bonds 
inclined  to  the  c-axis  than  are  parallel  to  it.  The 
results  of  this  averaging  are  given  in  Table  9,  along 
with  values  for  the  ideal  ice  structure  Ih,  but  we 
must  remember  that  the  defects  have  been  mod¬ 
eled  not  in  bulk  ice  but  in  a  cluster  of  only  eight 
molecules. 

The  first  thing  that  attracts  our  attention  is  that 
the  length  of  the  bond  carrying  the  defect  has 
increased.  This  increase  is  larger  for  L-defects, 
which  is  in  qualitative  agreement  with  the  larger 
activation  volume  for  L-defects  found  experimen¬ 
tally  (Table  7).  Plummer  discovered  a  large  in¬ 
crease  in  the  H-H  distance  in  a  D-defect.  One  of  the 
protons  actually  moves  off  the  hydrogen  bond, 
taking  an  interstitial  position,  and  the  remaining 
hydrogen  bond  becomes  bent  in  the  opposite  di¬ 
rection  by  approximately  20°. 

It  is  quite  remarkable  that,  in  spite  of  the  small 
size  of  simulated  clusters,  the  energy  of  pairs  of  L- 


Table  9.  Average  parameters  of  Bjerrum  defects 
in  ice. 


Parameter 

Defect 

Perfect 

structure 

D 

L 

0-0  distance  (A) 

2.865 

3.17 

2.76 

H-H  distance  (A) 

2.12 

- 

0.76 

Creation  energy  (eV) 

0.312 

0.336 

- 
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a.  Ionic  defect  in  one-dimensional  chain  of  hydrogen  bonds. 


Figure  11.  Soliton  model  of  OH  ion. 


and  D-defects  according  to  Plummer  is  about  0.65 
eV,  which  is  very  close  to  the  value  0.68  eV  found 
from  experiments. 

We  can  make  one  more  interesting  conclusion 
from  Plummer's  calculations  by  observing  that 
energies  of  defects  in  "eclipsed"  and  "staggered" 
configurations  are  different.  Because  in  diffusion 
and  electrical  conductivity  processes,  Bjerrum 
defects  must  flip  from  one  configuration  to  an¬ 


other,  the  difference  in  corresponding  energies 
may  contribute  to  the  observed  activation  energy 
of  defect  migration .  For  L-defects  this  difference  is 
only  0.026  eV  (i.e.,  on  the  order  oIIcrT),  whereas  for 
D-defects  it  is  equal  to  0.14  eV.  Is  this  not  the  reason 
for  lower  mobility  of  D-defects  compared  to  L- 
defects?  Calculations  of  the  type  performed  by 
Plummer  are  quite  complicated  and  have  not  yet 
been  made  for  ions. 


+<Po 


-<Po 


a.  Defect  in  one-dimensional  chain. 


Figure  12.  Soliton  model  of  Bjerrum  defect. 
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The  second  type  of  model  of  the  atomic  struc¬ 
ture  of  protonic  defects  is  the  so-called  soliton 
model,  in  which  the  defect  is  spatially  delocalized 
over  several  lattice  sites.  Both  electrostatic  and 
elastic  energies  of  the  defect  are  lowered  by  this. 
Originally  suggested  for  the  description  of  ion 
transfer  along  hydrogen  bond  chains  by 
Antonchenko  et  al.  (1983)  and  for  Bjerrum  defects 
by  Sergienko  (1986),  this  model  has  been  further 
developed  in  many  papers;  it  has  numerous  advo¬ 
cates  and  no  less  numerous  opponents.  A  detailed 
theoretical  description  of  this  model  can  be  found 
in  a  recent  book  by  Davydov  (1991).  A  simplified 
diagram  of  an  OH-  "soliton  ion"  in  a  linear  chain 
of  hydrogen  bonds  is  shown  in  Figure  11.  Notice 
that  all  protons  to  the  left  of  the  center  are  dis¬ 
placed  to  the  left  end  of  their  bonds  and  all  protons 
on  the  right  are  displaced  to  the  right.  During  the 
migration  of  such  an  ion  to  the  right,  the  protons 
on  the  bonds  are  transferred  to  the  left.  For  a  soliton 
to  arise  and  remain  stable,  at  least  two  conditions 
must  be  satisfied:  1)  the  protons  must  behave 
cooperatively,  i.e.,  they  must  strongly  interact  with 
each  other,  and  2)  a  strong  nonlinear  interaction  of 
the  proton  and  oxygen  sub-lattices  is  required. 

The  soliton  simulating  an  L-defect  (the  migra¬ 
tion  of  the  defect  from  left  to  right  transfers  the 
protons  from  right-side  positions  on  the  bonds  to 
left-side  positions)  is  shown  in  Figure  12.  Using 
the  same  schematic,  the  reader  can  easily  construct 
an  H30+  ion  and  a  D-defect.  Since  the  transfer  of 
soliton  defects  results  in  exactly  the  same  changes 
in  the  positions  of  protons  on  the  bonds  as  in  the 
case  of  transfer  of  their  "point"  analogs,  the  formal 
description  of  their  effects  in  the  ice,  effective 
charges,  contributions  to  ,  etc.,  will  be  exactly  the 
same. 

In  spite  of  a  certain  elegance  of  the  soliton 
model,  serious  arguments  have  appeared  against 
its  applicability  to  ice.  The  first  argument  pre¬ 
sented  by  Nagle  (1992)  is  that  solitons,  requiring 
for  their  motion  a  coordinated  displacement  of 
protons  arranged  into  some  ordered  positions, 
will  just  get  "stuck"  in  a  disordered  ice  protonic 
structure.  A  second  objection  is  that  the  soliton  is 
defined  and  has  to  move  along  a  particular  track, 
whereas  the  bonds  in  ice  form  a  branched  net¬ 
work.  The  relaxed  state  of  a  static  defect  may  in 
reality  be  somewhat  delocalized,  but  the  theory  of 
solitons  is  essentially  dynamic. 

A  third  objection  to  the  soliton  model  in  or¬ 
dered  water  molecule  chains  arises  from  precise 
ab  initio  calculations  performed  by  Godzik  (1990) 
and  Scheiner  (1991).  These  calculations  show  that 


in  ordered  chains  of  H2O  molecules  the  coopera¬ 
tive  motion  of  several  protons  involves  a  signifi¬ 
cantly  higher  energy  barrier  than  that  for  succes¬ 
sive  jumps  of  one  proton  after  another. 

We  cannot  leave  the  structure  of  protonic  de¬ 
fects  without  mentioning  the  possibility  that  these 
defects  could  be  coupled  with  vacancies  (Kroger 
1974).  For  example  a  D-defect  at  a  vacancy  is 
merely  a  vacant  site  with  three  rather  than  the 
normal  two  hydrogen  bonds  pointing  into  it,  and 
its  D  character  will  move  with  it  as  it  moves 
through  the  lattice  (see  Fig.  17).  As  vacancies  may 
be  highly  mobile,  we  must  allow  for  the  possibility 
that  the  majority  of  Bjerrum  defects  could  exist 
and  move  in  this  way.  There  is,  however,  no  ex¬ 
perimental  evidence  on  this  issue. 

IMPURITIES 

Although  liquid  water  is  a  very  good  solvent  for 
many  substances,  its  solid  phase — ice — dissolves 
almost  nothing.  Most  of  the  impurities  dissolved 
in  water  either  get  rejected  as  the  water  freezes  or 
precipitate  as  a  second  phase  within  the  ice.  There 
are  some  exceptions,  such  as  a  few  acids  (HF,  HC1), 
ammonia  (NH3),  some  alkalis  (KOH,  NaOH)  and 
their  derivatives  like  NH4F.  These  substances  can 
be  incorporated  properly  into  the  ice  structure 
where  they  drastically  change  the  protonic  carrier 
concentrations  in  the  ice.  To  understand  how  these 
impurities  affect  the  electrical  properties  of  ice,  let 
us  consider  the  probable  patterns  of  their  incorpo¬ 
ration  into  the  ice  lattice  (see  Fig.  13-15).  Note  that 
although  there  now  exist  a  number  of  experimen¬ 
tal  techniques  for  investigating  crystal  structures 
on  an  atomic  scale  (high  resolution  transmission 
electron  microscopy,  scanning  tunneling  micros¬ 
copy,  atomic  force  microscopy,  etc.),  such  tech¬ 
niques  have  not  yet  been  successfully  applied  to 
ice.  Patterns  such  as  those  represented  in  Figures 
13-15  turn  out  to  be  just  hypotheses  based  upon 
ionic  radii,  energies  and  lengths  of  atomic  bonds, 
etc.  The  ionic  radius  of  F-,  for  example,  is  rather 
close  to  that  of  O2-,  and  it  is  natural  to  assume  that 
fluorine  substitutes  for  oxygen  successfully  in  an 
ice  lattice.  As  can  be  seen  from  Figure  13,  since  the 
UF  molecule  can  provide  only  one  proton  for  the 
four  nearest  hydrogen  bonds,  one  bond  appears  to 
be  lacking  any  protons  at  all;  i.e.,  substitution  of 
the  water  molecule  by  the  HF  molecule  results  in 
the  appearance  of  an  L-defect  associated  with  the 
impurity  molecule.  The  binding  energy  of  such  an 
L-defect  with  the  HF  molecule  is  not  large  (<  10-2 
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eV  [Jaccard  1959,  Camplin  et  al.  1978]).  Since  the 
number  of  possible  arrangements  of  the  L-defect 
in  an  ice  lattice  is  enormous,  practically  all  L- 
defects  leave  their  parent  molecules. 

Moreover,  it  turns  out  that  the  single  proton  in 
the  HF  molecule  is  not  bound  firmly  and  can  jump 
along  a  hydrogen  bond  to  another  water  molecule, 
turning  that  water  molecule  into  an  H3O  ion.  It  is 
important  to  emphasize  the  chemistry  of  HF.  To  a 
physicist  the  natural  thing  would  be  for  an  extra 


proton  to  attach  to  HF,  forming  an  OH-  ion,  but 
chemically  HF  is  a  weak  acid  in  water  and  prefers 
to  lose  a  proton. 

At  the  same  time  as  HF  doping  of  ice  increases 
the  L-defect  and  HaO+  ion  concentrations,  it  will 
also  lower  the  concentrations  of  D-defects  and 
OH-  ions.  This  happens  because  the  equilibrium 
concentrations  of  D-defects  and  OH-  ions  in  pure 
ice  are  determined  by  a  balance  between  the  rate  of 
their  thermoactivated  production  and  the  rate  of 


a.  Ice  lattice  with  incorporated  HF  molecule. 


b.  Same  after  HgO+  ion  and  L-defect  were  released  by  means  of  sequential  hops  of  protons 
(along  the  hydrogen  bonds  to  release  HjO+  ion  and  between  the  bonds  to  release  L-defect). 

Figure  13.  Hypothesized  pattern  of  incorporation  ofHF  impurity. 
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recombination:  E>-  with  L-defects  and  H3O+  with 
OH-  ions.  HF  doping  increases  L-defect  and  H3O* 
ion  concentrations,  so  that  the  rate  of  recombina¬ 
tion  increases,  which  results  in  decreasing  D-de- 
fect  and  OH-  ion  concentrations.  This  can  be  stated 
in  terms  of  the  law  of  mass  action,  which  is  that  the 
products  of  concentrations  M|Xn2  and  M3  x  M4  must 
be  constant  in  thermal  equilibrium,  irrespective  of 
doping 

»1  x  m2  =  -4  (Mujofcxp  J-  ,  «i,2  <<  «it20  ^69) 


”3  X  ”4  =  (MH2ofexP  /  ”3,4  «  ”HzO  :  (70) 

Another  molecule  that  can  successfully  be  in¬ 
corporated  into  the  ice  lattice,  substituting  for 
water  molecules,  is  ammonia,  NH3  (Fig.  14).  Since 
NH3  donates  three  protons  (instead  of  two  as  in 
the  case  of  H2G)  for  four  hydrogen  bonds,  we 
observe  two  protons  at  one  bond  at  the  same  time, 
i.e.?  a  D-defect.  As  a  consequence  of  the 
thermoactivation  process,  this  defect  can  be  re¬ 
leased  and  fake  part  in  the  conductivity  process. 


a.  NH3  molecule  incorporated  into  ice  structure. 


b-  Same  after  release  of  OH  ions  and  D-defect  from  the  NH3  molecule. 


figure  14.  Hypothesized  pattern  of  incorporation  ofNH 3  impurity. 
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Besides,  the  NH3  molecule  is  a  base  and  can  accept 
a  fourth  proton,  forming  an  OH-  ion  at  a  neighbor¬ 
ing  water  molecule.  Activation  energies  for  the 
release  of  D-defects  and  OH-  ions  from  the  ammo¬ 
nia  molecule  in  ice  are  comparatively  large,  so 
that,  although  ammonia  doping  of  ice  increases  D- 
defect  and  OH-  ion  concentrations  (decreasing  L- 
defect  and  H3O  ion  concentrations  in  accordance 
with  eq  69-70),  its  influence  in  terms  of  conductiv¬ 


ity  is  not  as  efficient  as  that  of  HF.  It  is  also 
important  that  NH3  doping  increases  the  concen¬ 
tration  of  less  mobile  carriers  (D-defects  and  OH- 
ions);  see  Table  6. 

Lastly,  the  third,  very  important  group  of  im¬ 
purities  that  considerably  affect  the  conductivity 
of  ice  are  the  alkali  hydroxides,  among  which  the 
most  significant  is  KOH.  The  most  plausible  mecha¬ 
nism  of  incorporation  of  KOH  into  the  ice  struc- 


a.  Hoiv  KOH  dissolved  in  ice  is  thought  to  introduce  defects.  The  K+  ion  is  interstitial,  the 
OH~  ion  substitutes  for  a  water  molecule  leaving  one  bond  without  any  proton — an  L-defect. 


b.  Same  after  L-defect  and  OH  ion  have  moved  away  from  the  K+  ion. 
Figure  15.  Hypothesized  pattern  of  incorporation  of  KOH  impurity. 
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Figure  16.  Absorption  spectra  of  solvated  electrons  in 
water  and  ice  (after  Nilsson  et  al.  1972). 


ture  is  shown  in  Figure  15.  The  potassium  atom  is 
accommodated  in  interstices,  while  the  hydroxyl 
group  OH  is  built  into  the  lattice.  As  can  be  seen 
from  Figure  15,  this  results  in  the  appearance  of  L- 
defects  and  OH~  ions  linked  to  KOH  impurity. 
Thermoactivation  processes  allow  them  to  be  re¬ 
leased  into  the  ice  bulk,  increasing  L-defect  and' 
OH-  ion  concentrations  and  suppressing  D-de- 
fects  and  H30+  ions.  Because  of  the  extremely 
small  binding  energies  of  L-defects  and  OH-  ions 
with  K+  ions  in  ice,  these  impurities  become  effi¬ 
cient  donors  of  L-defects  and  OH-  ions  even  at 
very  low  temperatures  (Zaretskii  et  al.  1988). 

The  subject  of  the  concentrations  of  charge  car¬ 
riers  in  doped  ice  has  been  treated  many  times  in 
the  scientific  literature  (see  Jaccard  1959,  Kroger 
1974,  Camplin  et  al.  1978),  and  is  described  more 
fully  in  the  report  on  the  Electrical  Properties  of  Ice 
(Petrenko  1993a).  For  the  fullest  detail,  based  on 
the  law  of  mass  action,  the  reader  is  referred  to 
Kroger  (1974). 

At  present  there  are  hundreds  of  papers  avail¬ 
able  in  the  literature  in  which  the  electrical  proper¬ 
ties  of  doped  ice  were  studied.  The  results  of  these 
studies  were  discussed  in  detail  in  the  review  by 
Petrenko  (1993a)  and  are  presented  in  part  in 
Tables  4-6.  Let  us  refer  to  the  following  works  as 


examples:  Takei  and  Maeno  (1984, 1987)  on  HC1 
doped  ice;  Gross  (1975)  and  Gross  et  al.  (1978)  on 
HC1,  HF,  NaCl,  KF,  NH4F,  NH4OH,  NH4CI, 
NH4HCO3  doped  ice;  KOH  doped  ice  was  studied 
by  Zaretskii  et  al.  (1988)  and  by  Howe  and 
Whitworth  (1989). 

On  the  whole,  the  effect  of  impurities  dissolved 
in  ice  on  its  electrical  properties  is  rather  well 
explained  on  the  basis  of  the  currently  accepted 
model  of  electrical  conductivity  of  ice,  provided 
that  the  concentrations  of  all  four  types  of  charge 
carriers  are  calculated  following  the  method  sug¬ 
gested  by  Kroger  (1974). 

The  impurity  diffusion  coefficients  in  ice  were 
measured  in  a  number  of  works.  Haltenorth  and 
Klinger  (1969, 1977)  investigated  diffusion  of  HF 
in  monocrystalline  ice  parallel  and  perpendicular 
to  the  c-axis.  They  found  that  at  T  =  -10°C  the 
diffusion  coefficient  of  HF  in  ice  is  1.08  x  10-11  m2 
s-1  and  obeys  the  Arrhenius  law  with  an  activation 
energy  of  0.2  eV.  The  value  of  the  diffusion  coeffi¬ 
cient  along  the  c-axis  was  20%  lower  than  the 
corresponding  value  normal  to  it.  In  polycrystal¬ 
line  ice  the  diffusion  coefficient  was  25%  higher 
than  in  monocrystals  owing  to  grain  boundary 
diffusion.  It  should  be  noted  that  the  diffusion 
coefficient  for  HF  is  much  larger  than  the  self¬ 
diffusion  coefficients  of  H2O  given  in  Table  2, 
which  suggests  that  HF  molecules  move  through 
interstitial  positions  more  easily  than  those  of 

h2o. 

Diffusion  of  the  neutral  gases  He  and  Ne  was 
studied  by  Haas  et  al.  (1973).  Diffusion  coefficients 
of  these  gases  in  ice  appeared  to  be  extremely  high, 
indicating  a  "transparency"  of  ice  for  these  gases. 
It  is  interesting  that  the  diffusion  coefficients  along 
the  c-axis  (D|  |c)  were  six  times  (He)  and  eight 
times  (Ne)  higher  than  in  the  perpendicular  direc¬ 
tion.  In  the  temperature  range  173-253  K  for  He 
and  213-263  K  for  Ne,  Haas  and  co-authors  found 
that  the  following  relationships  hold 

D 1  ic(He)  =  3.4  x  10-7  exp  m2  s-1  (71) 

Dlc(He)  =  1.1  x  10-7  exp  j  m2  s_1  (72) 

D  |  ic(Ne)  =  1.2  x  10-6  exp  f-  0-25  eV\  m2  s-l  (73) 

ELECTRONIC  DEFECTS 

In  this  section  we  will  consider  point  defects 
that  arise  from  an  excitation  of  the  electronic  sub¬ 
system  of  ice:  hydrated  electrons,  radicals  and 
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Figure  17.  Complex  (vacancy  +  D-defect)  formation  in  ice  (after  de  Haas  et  ah  1983). 
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color  centers.  Since  ice  has  a  very  wide  forbidden 
band  with  Eg  =  8.5  eV  (Ryzhkin  1992),  the  above 
excitations  usually  arise  when  ice  is  either  illumi¬ 
nated  by  hard  ultraviolet  radiation  or  is  excited  by 
beams  of  charged  particles,  y-rays  and  X-rays. 

Conducting  electrons  in  ice  have  been  studied 
in  many  works  (e.g.,  Verbeme  et  al.  1978;  Warman 
et  al.  1980, 1982, 1983;  de  Haas  et  al.  1983;  Kunst  et 
al.  1983).  Conducting  electrons  in  these  studies 
were  excited  either  by  a  short  X-ray  pulse  of  about 
0.5  ns,  produced  by  irradiation  of  a  Pt  target  with 
3-MeV  electrons  from  a  Van  de  Graaf  accelerator, 
or  directly  by  the  3-MeV  electron  beam.  The  life¬ 
time  of  electrons  in  the  conduction  band  is  ex¬ 
tremely  small;  at  T  =  -20°C  it  is  several  nanosec¬ 
onds  and  at  T= -60°C  several  tens  of  nanoseconds 
(de  Haas  et  al.  1983).  The  main  processes  limiting 
the  lifetime  of  electrons  at  these  temperatures  are 
trapping  and  solvation  (hydration).  In  spite  of 
small  lifetimes,  the  analysis  of  transient  electron 
currents  in  strong  electric  fields  showed  their 
mobility  to  be  (ie  =  (25  ±  5)  x  ICHm2  V”1  s-1  (Ver¬ 
beme  et  al.  1978,  Warman  et  al.  1980, 1982).  This 
value,  as  the  reader  must  have  noticed  before,  ex¬ 
ceeds  the  mobility  of  most  mobile  protonic  charge 
carriers  in  ice  by  more  than  four  orders  of  magni¬ 
tude  (see  Table  5).  The  mobility  of  electrons  in  ice 
is  almost  independent  of  temperature  from  -60 
down  to  -120°C. 

All  conducting  electrons  move  very  rapidly 
into  the  so-called  solvated  state.  In  the  literature 
these  electrons  are  often  referred  to  as  "solvated" 
or  "hydrated"  electrons;  these  terms  are  equiva¬ 
lent.  Such  electrons  manifest  themselves  through 
a  specific  absorption  spectrum  with  a  maximum  at 
approximately  6800  A  and  through  a  distinct  Elec¬ 
tron  Spin  Resonance  (ESR)  line.  Initially,  solvated 
electrons  were  discovered  in  water  and  later  in  ice. 
A  great  number  of  publications  are  dedicated  to 
solvated  electrons  and  yet  many  of  their  proper¬ 
ties  remain  unclear.  Below,  we  shall  address  ques¬ 
tions  about  the  structure,  lifetime,  optical  proper¬ 
ties,  recombination  and  mobility  of  solvated  elec¬ 
trons  in  ice. 

The  major  reason  for  electrons  first  to  become 
localized  and  then  to  move  into  the  solvated  state 
is  considered  to  be  the  polarization  of  the  sur¬ 
rounding  medium  by  the  electron.  Owing  to  a 
large  dipole  moment,  water  molecules  tend  to 
orient  towards  the  electron.  The  electrostatic  en¬ 
ergy  of  the  state  in  which  the  electron  is  "coated" 
by  water  molecules  oriented  towards  it  would  be 
obviously  less  than  the  initial  nonpolarized  state. 
In  water,  where  each  molecule  can  be  reoriented 


extremely  fast,  such  polarization  would  seem  to 
occur  in  a  different  fashion  than  in  ice  in  which  the 
mutual  molecular  orientation  follows  the  ice  rules. 
Nevertheless,  the  optical  absorption  spectrum  of 
the  final,  fully  relaxed  state  of  the  trapped  electron 
in  ice  resembles  very  closely  that  of  the  solvated 
electron  in  water,  with  only  a  small  shift  in  the 
position  of  the  absorption  maximum  at  1.8  eV 
being  found  upon  melting  (Nilsson  et  al.  1972, 
Shubin  et  al.  1966,  Taub  and  Eiben  1968). 

Several  of  such  spectra  in  ice  and  one  in  water 
are  shown  in  Figure  16.  The  observation  of  similar 
absorption  spectra  in  water-containing  media  can 
be  considered  as  evidence  for  the  presence  of 
solvated  electrons.  This  similarity  of  the  absorp¬ 
tion  spectra  of  solvated  electrons  in  water  and  ice 
strongly  suggests  the  same  basic  microscopic  struc¬ 
ture  for  the  solvated  state  in  both  media.  Present 
opinion,  based  on  theoretical  considerations  and 
ESR  studies  of  electrons  in  low-temperature  aque¬ 
ous  glasses  (Kevan  1981),  tends  to  favor  a  cavity 
model  for  the  solvated  electron,  i.e.,  a  model  in 
which  the  electron  occupies  a  void  surrounded  by 
oriented  water  molecules.  This  being  the  case,  the 
initial  localization  process  in  ice  must  therefore 
involve  trapping  at  a  vacancy  in  the  lattice. 

Figure  17,  adopted  from  the  paper  by  de  Haas  et 
al.  (1983),  shows  schematically  one  of  the  possible 
structures  of  the  solvated  electron  in  ice.  Such  a 
defect  is  a  combination  of  a  water  molecule  va¬ 
cancy,  a  D-defect  and  an  electron  (not  shown).  In 
this  type  of  a  defect,  the  electrostatic  energy  is 
lowered  by  attraction  between  the  electron  and 
three  protons.  A  migration  of  this  defect  must  be 
accompanied  by  cooperative  drift  of  the  vacancy 
and  the  D-defect.  Their  mobilities  could  have  pro¬ 
vided  some  estimation  of  the  solvated  electron's 
mobility  in  ice.  Unf  ortunately,  we  could  not  find  in 
the  literature  any  reliable  data  on  diffusion  coeffi¬ 
cients  of  either  vacancies  or  D-defects  or  solvated 
electrons  in  ice.  The  mobility  of  solvated  electrons 
in  water  was  measured  to  be  1.84X10-7  m2  V-1  s-1 
(Schmidt  and  Buck  1966).  This  value  can  be  con¬ 
sidered  as  an  upper  limit  for  mobility  of  solvated 
electrons  in  ice. 

Both  in  water  and  ice,  solvated  electrons  are 
annihilated  by  the  recombination  reaction  (Muto 
et  al.  1992) 

e-  +  H3O+  ->  H20  +  H*  •  (74) 

In  addition,  they  can  react  with  impurities  in 
the  bulk  and  at  the  surface.  All  these  factors  reduce 
their  lifetimes,  which  nonetheless  can  be  very  large 
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at  low  temperatures  (about  80  minutes  at  4.2  K 
[Muto  et  al.  1992]). 

In  eq  74  we  have  used  the  symbol "  * "  to  denote 
a  "radical,"  i.e.,  an  atom  or  a  molecule  that  has  an 
unpaired  electron  spin.  Owing  to  the  presence  of 
such  an  electron,  the  radicals  are  extremely  chemi¬ 
cally  active  and  consequently  unstable.  In  eq  74  we 
have  a  hydrogen  radical  or  simply  a  hydrogen 
atom  H*;  H*  is  not  the  only  radical  produced 
during  ionization  of  water  molecules  in  ice.  We  list 
below  some  most  frequently  encountered  reac¬ 
tions  in  which  radicals  consisting  of  hydrogen  and 
oxygen  atoms  are  created  and  annihilated: 


H2O  +  hv  — >  H20+  +  e~ 

(75) 

H20+  +  H20  ->  H30+  +  OH* 

(76) 

H3O+  +  e-  -»  H20  +  H* 

(77) 

OH*  +  OH*  -»  H202 

(78) 

H*  +  H*  -»H2 

(79) 

H*  +  OH*  H20 

(80) 

OH*  +  H20  s  ±5  H30+  +  O- 

(81) 

H202  +  H*  ->  H20  +  OH* 

(82) 

H202  +  OH*  ->  H02*  +  H20  . 

(83) 

The  study  of  the  properties  of  radicals  belongs  to 
the  extensive  domain  of  radiation  chemistry  and 
falls  outside  the  scope  of  this  report.  The  reader 
interested  in  the  properties  of  radicals  in  ice  can 
find  a  detailed  list  of  references  on  this  question  in 
the  papers  by  Bednarek  and  Plonka  (1987a, b)  and 
also  in  the  already  cited  references  concerning  free 
(conduction)  and  solvated  electrons. 

Let  us  indicate  here  that,  because  of  the  frac¬ 
tional  charge  of  ions  in  ice  (±  0.62e),  some  of  the 
reactions  (eq  75-83)  (written  for  water,  where  ions 
have  an  integer  charge)  must  be  altered  and 
completed  to  fulfill  the  charge  conservation  law. 
For  example,  the  charge  of  particles  in  the  left- 
hand  side  of  eq  77  is  -0.38e,  and  that  of  the  right- 
hand  side  is  zero.  The  charge  conservation  law  will 
hold  if  we  add  to  the  left-hand  side  a  D-defect.  In 
this  case  the  electron  is  a  solvated  electron  trapped 
by  the  complex:  a  water  molecule  vacancy  +  D- 
defect.  Thus,  now  the  reaction  (eq  74)  can  be 
written  as 

H30+  +  (e-  +  VH2O  +  D)  -»  H2O  +  H*  +  VH2O  (84) 


For  the  charge  to  be  conserved  in  the  above  reac¬ 
tion  (eq  77),  we  have  either  to  add  a  D-defect  to  the 
right-hand  side  or  to  assume  that  the  radical  OH* 
has  a  charge  of  +0.38e.  The  possibility  for  the 
radicals  in  ice  to  be  electrically  charged  was  first 
suggested  by  Chesnakov  et  al.  (1987).  The  reasons 
for  the  radicals  embedded  in  the  ice  structure  to 
have  an  electric  charge  are  the  same  that  lead  to  the 
fractional  charges  of  ions  in  ice.  There  is  no  experi¬ 
mental  evidence  available  that  supports  or  refutes 
the  presence  of  electrical  charge  at  radicals  in  ice. 
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